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Abstract 

We consider random Schrodinger equations on M. d for d > 3 with a homogeneous 
Anderson-Poisson type random potential. Denote by A the coupling constant and the 
solution with initial data po- The space and time variables scale as x r^i A -2—K / 2 , t ~ 

A -2 ~ k with 0 < k < Ko(d). We prove that, in the limit A —> 0, the expectation of the 
Wigner distribution of pj converges weakly to the solution of a heat equation in the 
space variable x for arbitrary L 2 initial data. 

The proof is based on analyzing the phase cancellations of multiple scatterings on 
the random potential by expanding the propagator into a sum of Feynman graphs. In 
this paper we consider the non-recollision graphs and prove that the amplitude of the 
non-ladder diagrams is smaller than their “naive size” by an extra A c factor per non- 
(anti)ladder vertex for some c > 0. This is the first rigorous result showing that the 
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improvement over the naive estimates on the Feynman graphs grows as a power of the 
small parameter with the exponent depending linearly on the number of vertices. This 
estimate allows us to prove the convergence of the perturbation series. 

AMS 2000 Subject Classification: 60J65, 81T18, 82C10, 82C44 

1 Introduction 

The fundamental equations governing the basic laws of physics, the Newton and the 
Schrodinger equations, are time reversible and have no dissipation. It is remarkable 
that dissipation is nevertheless ubiquitous in nature, so that almost all macroscopic 
phenomenological equations are dissipative. The oldest such example is perhaps the 
equation of heat conductance found by Fourier. 

On a microscopic level, Brown observed almost two centuries ago that the motion 
of a pollen suspended in water was erratic [5j. This led to the kinetic explanation 
by Einstein in 1905 [14] that such a motion was created by the constant “kicks” on 
the relatively heavy pollen by the light water molecules. It should be noted that at 
that time even the atomic-molecular structure of matter was not universally accepted. 
Einstein’s theory was strongly supported by Boltzmann’s kinetic theory, which, however, 
was phenomenological and seriously debated at the time. Finally in 1908 Perrin [40 | 
experimentally verified Einstein’s theory and used it, among others, to give a precise 
estimate on the Avogadro number. These experiments gave the strongest evidence for 
atoms and molecules at that time. 

In Einstein’s kinetic theory both the heavy particle (the pollen) and the light par¬ 
ticles (the water molecules) obey Newton’s law. Therefore, Einstein’s kinetic theory in 
fact postulated the emergence of the Brownian motion from a classical non-dissipative 
reversible dynamics. The key difficulty of a mathematically rigorous derivation of Brown¬ 
ian motion from reversible dynamics is similar to the justification of Boltzmann’s molec¬ 
ular chaos assumption (StoBzahlansatz); the dissipative character emerges only in a 
scaling limit, as the number of degrees of freedom goes to infinity. 

The first mathematical definition of the Brownian motion was given in 1923 by 
Wiener, who constructed the Brownian motion as a scaling limit of random walks. This 
construction was built upon a stochastic microscopic dynamics which by itself is dissipa¬ 
tive. The derivation of the Brownian motion from a time reversible Hamiltonian system, 
however, was not seriously considered until more than half a century later. Kesten- 
Papanicolaou |33] in 1978 proved that the velocity distribution of a particle moving in 
an environment consisting of random scatterers (i.e., Lorenz gas with random scatterers) 
converges to a Brownian motion in a weak coupling limit for d > 3. In this model the 
bath of light particles is replaced with random static impurities. The same result was 
obtained later in d = 2 dimensions by Durr, Goldstein and Lebowitz m- In a recent 
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work [26], Komorowski and Ryzhik have controlled the same evolution on a longer time 
scale and proved the convergence to Brownian motion of the position process as well. 
Bunimovich and Sinai [9] proved the convergence of the periodic Lorenz gas with a hard 
core interaction to a Brownian motion in 1980. In this model the only source of random¬ 
ness is the distribution of the initial condition. Finally, Durr, Goldstein and Lebowitz 
m proved that the velocity process of a heavy particle in a light ideal gas, which is a 
model with a dynamical environment, converges to the Ornstein-Uhlenbeck process. 

Wiener’s construction of Brownian motion is based on a random walk. The ran¬ 
dom walk could easily be replaced by the Markovian process generated by a linear 
Boltzmann equation. The linear Boltzmann equation was rigorously derived from the 
classical Lorenz gas by Gallavotti [29], Spohn |36j and Boldrighini, Bunimovich and 
Sinai 0. (The nonlinear Boltzmann equation was derived by Lanford [[37] for short 
time.) Although Brownian motion was discovered and theorized in the context of clas¬ 
sical dynamics, we shall prove that it also describes the motion of a quantum particle in 
a random environment. 

The random Schrodinger equation, or the quantum Lorentz model, is given by the 
evolution equation: 

id t ^ t {x) = Hi/j t (x), H = + XVu(x) , (1.1) 

where A is the coupling constant and is the random potential. The first scale with a 
non-trivial limiting dynamics is the weak coupling limit, A —> 0, where the space, time 
and the coupling constant are subject to the kinetic scaling: 

t —> te _1 , x —> xe -1 , A = \fe . (1.2) 

Under this limit, the appropriately rescaled Wigner distribution (see (12.111) 1 of the solu¬ 
tion to the Schrodinger evolution m converges weakly to a linear Boltzmann equation. 
This was first established by Spohn [45] for a random potential with Gaussian distri¬ 
bution and small macroscopic time. This method was extended to study higher order 
correlations in [28]. A different method (applicable to the lattice setting and general 
random potential, see remarks later on) was developed in [T6] where the short time 
restriction was removed. This method was also extended to the phonon case m and 
Lukkarinen and Spohn [39] have employed a similar technique for studying the energy 
transport in a harmonic crystal with weakly perturbed random masses. 

Since the long time limit of a Boltzmann equation is a heat equation, we shall take a 
time scale longer (see ([2.171) 1 than in the kinetic scaling limit (|1.2I) . Our aim is to prove 
that the limiting dynamics of the Schrodinger evolution in a random potential under 
this scaling is governed by a heat equation. This requires to control the Schrodinger 
dynamics up to a time scale A~ 2 ~ K , k > 0. Quantum correlations that are small on the 
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kinetic scale and are neglected in the first limit may contribute on the longer time scale. 
The derivation of the heat equation is thus much more difficult than first deriving the 
Boltzmann equation from Schrodinger dynamics on the kinetic scale and then showing 
that Boltzmann equation converges to a diffusive equation under a different limiting 
procedure. Notice that the limit in our approach is a long time scaling limit which 
involves no semiclassical limit. 

The approach of this paper also applies to lattice models and yields a derivation of 
Brownian motion from the Anderson model PH 18] and in fact, we present our main 
technical steps in a unified framework. The dynamics of the Anderson model was postu¬ 
lated by Anderson to be localized for large coupling constant A and extended for small 
coupling constant (away from the band edges and in dimension d > 3). The localization 
conjecture was first established rigorously by Goldsheid, Molchanov and Pastur [32] in 
one dimension, by Frohlich-Spencer [25], and later by Aizenman-Molchanov [T] in several 
dimensions, and many other works have since contributed to this field. The progress for 
the extended state conjecture, however, has been limited. It was proved by Klein [MJ 
that all eigenfuctions are extended on the Bethe lattice (see also [2, 24J). In Euclidean 
space, Schlag, Shubin and Wolff [53] proved that the eigenfunctions cannot be localized 
in a region smaller than A~ 2+5 for some 6 > 0 in d = 2. Chen HQ], partly based on 
the method [16], extended this result to all dimensions d> 2 and 5 = 0 (with logarith¬ 
mic corrections). Extended states for Schrodinger equation with a sufficiently decaying 
random potential were proven by Rodnianski and Schlag [SI] and Bourgain [8j (see also 

HD)- 

In summary, all known results for the Anderson model (or its modifications) in Eu¬ 
clidean space are in regions where the dynamics have either no effective collision or there 
are typically only finitely many of them. Under the diffusive scaling of this paper, see 
(12.171) . the number of effective scatterings is a negative fractional power of the scaling 
parameter. In particular, it goes to infinity in the scaling limit, as it should be the case 
if we aim to obtain a Brownian motion. As in [10] . one may derive from our dynamical 
result that the eigenfunctions cannot be localized in a region smaller than \~ 2 ~ K / 2 and 
dimension d > 3. 

Acknowledgement. The authors are grateful for the financial support and hospitality 
of the Erwin Schrodiger Institut, Vienna, Max Planck Institut, Leipzig, Stanford Uni¬ 
versity and Harvard University, where part of this work has been done. We thank Jani 
Lukkarinen for pointing out a simplification of Lemma 16.11 for the continuum model and 
we also thank the referee for his comments and suggestions. 
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2 Statement of the main result 


2.1 Notations 

Let 

H := -^-A + XV (2.1) 

2 m 

denote a random Schrodinger operator acting on L 2 (K d ), d > 3, with a random potential 
V = V UJ (x ) and a small positive coupling constant A. We shall choose the units so that 
h = 1, m = 1 so that fr 2 /(2m) = [2(27t) 2 ] _1 . The potential is given by 

K,(x) := [ B(x - y)dy UJ (y) , (2.2) 

JR d 

where B is a single site potential profile and y u is a Poisson point process on with 
homogeneous unit density and with independent, identically distributed random masses. 
More precisely, for almost all realizations u consists of a countable, locally finite collection 
of points, (y 7 (w) G : 7=1,2,...}, and random weights {u 7 ((u) G M : 7=1,2,...} 
such that the random measure is given by 

OO 

fiu = , ( 2 - 3 ) 

7=1 

where 5 y denotes the Dirac mass at y G M d . The Poisson process {y 7 (cu)} is independent 
of the weights {u 7 (cu)}. The weights are real i.i.d. random variables with distribution 
P, and with moments rrik := satisfying 

?n 2 = 1, rri 2 d < oo, mi = m 3 = m 5 = 0. (2.4) 

The expectation with respect to the random process {y 7 ,u 7 } is denoted by E. 

For the single-site potential, we assume that B is a spherically symmetric Schwarz 
function with 0 in the support of its Fourier transform, i.e. 

0 G supp (B) . (2.5) 

More precisely, we introduce the norm 

ll/IL,» := £ ||<i)”a"/(i)||oo 

|a|<n 

with (x) := (1 + x 2 ) 1 / 2 (here a is a multiindex) and we assume 

||-B||fc,fc < Ck Vfe G N . (2.6) 
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Actually, it is sufficient to assume (12.61) for all k < ko(n). 

We note that the operator H w is not bounded from below due to the possible large 
concentration of Poisson points in some region. Nevertheless, H w is self-adjoint under 
very general conditions, see [35] . 

We introduce a few notational conventions. The letters x,y,z will always denote 
configuration space variables, while p , q, r, u, v, w will be reserved for d-dimensional mo¬ 
mentum variables. The norm without indices, || • ||, will always denote the standard 
L 2 (R d ) norm. The bracket denotes the standard scalar product on L 2 (R d ) and 
(■, ■} will denote the pairing between the Schwarz space and its dual on the phase space 
R d x R d . 

Integrals without explicit domains will always denote integration over R d with respect 
to the Lebesgue measure. For any / E L 2 (R d ), the Fourier transform is given by 

Tip) ■■= J e~ 2 ^ x f(x)dx , P eR d . (2.7) 

and the inverse Fourier transform is given by 

g(x) = J g(p)e 2mp - x dp , x E R d . 

For functions defined on the phase space, f(x,v), the Fourier transform will always be 
taken only in the space variable, i.e. 

M, v) ■■= J e~ 2 ^ x f(x, v)dx , (6^. 

The Fourier transform of the kinetic energy operator is given by 

[ - A/] (p) = [ - I • pAA/] (p) = , 

where 

<p)-=\p 2 ( 2 - 8 ) 

is the dispersion law. The velocity is given by ^Ve(p) = 4^ P- We note that the addi¬ 
tional h 2 /m = (2 tt)- 2 factor was erroneously omitted in the definition of the Hamiltonian 
in (1.2) of [T9]. 

Define the Wigner transform of a function if E L 2 (R d ) 

W,p(x,v) := J e 2mv ' v il;(x + ^)V>(® - ^)drj . (2.9) 
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( 2 . 10 ) 


The Fourier transform of W^(x, v) in the x variable is therefore 

W^{£,v) = $(v - |)^(u + |) . 

Define the rescaled Wigner distribution as 

W$(X,V) :=e~ d W^,v) . ( 2 . 11 ) 

Its Fourier transform in X is given by 

w$&v) = $(y-^)$(y + ^) . 

For any function h : —> C and energy value e > 0 we introduce the notation 

[h](e) ■= J h(v)S(e - e(v))dv := J ^ h(q) ^(q)\ ’ ( 2 - 12 ) 

where du{q) is the restriction of the Lebesgue measure onto the energy surface £ e : = 
{q : e(q) = e} that is the ball of radius \/2e. More explicitly, 

[h](e) := ( 2 e) 2 _1 f Zi(\/2e0)d0 . 

Jsd-i 


Clearly 


• /*oo 

h{v)dv = / [h\(e)de . 

Jo 


The normalization of the measure [-] e is given by 

[1]( e ) :=c d _ 1 (2e)|- 1 , 

where c^-i is the volume of the unit sphere S d ~ 1 . 


(2.13) 

(2.14) 


2.2 Main Theorem 

The weak coupling limit is defined by the following scaling: 

t = T/e , x = X/e, e = A 2 . (2-15) 

The Wigner distribution (X,V) converges weakly to a function Fq-(X, V) that 

satisfies the Boltzmann equation 

d T F T (X, V) + V ■ V x F r {X, V ) 

= 2vr J dU\B(U -V)\ 2 5{e{U) - e(V)) F T (X, U) - F r {X, V) . (2.16) 


7 



Note that the Boltzmann equation can be viewed as the generator of a Markovian semi¬ 
group on phase space. The proof of (I2.16P for continuum Gaussian model was given in 
[16] : the Z d lattice case with general i.i.d. random potential was considered in [TO]- The 
derivation of the Boltzmann equation for potential (12.21) follows from these two proofs 
in a straightforward way. 

In this paper we consider the long time scaling 

t = \- K (X~ 2 T), x = \- K/2 (\~ 2 X) = X/e, £ = \ k/2+2 (2.17) 

with some k > 0. This scaling corresponds to the long time limit of the Boltzmann 
equation with diffusive scaling. 

For some energy e > 0, let 

L ef{v) ■= J du a(u,v)[f(u) - f{v)\, e(v)=e, (2.18) 

be the generator of the momentum jump process v(t) on E e with collision kernel 

a(u, v ) := 2ir\B(u — v)\ 2 5(e(u) — e(v)) . (2.19) 

Notice that we use u, v etc. to denote the momenta; the corresponding velocities are 
u/(2 tt), v/(2n) etc. 

A well-known argument shows that B ^ 0 and the regularity of B guarantees the 
following properties. Some details will be given in Pi- 

Lemma 2.1 For each e > 0 the Markov process {u(f)}t>o with generator L e is uni¬ 
formly exponentially mixing. The unique invariant measure is the uniform distribution, 
[• ](e)/[l](e), on the energy surface £ e . 

Let 

D ij( e )'■= 77 TT 2 £ e [v W (t)u O) (0)]df , v = (u (1) ,... ,u (d) ), i,j = 1,2,. ..d, 

l Z7r J Jo 

be the velocity autocorrelation matrix, where £ e denotes the expectation with respect 
to this Markov process in equilibrium. By the spherical symmetry of B and e(U), the 
autocorrelation matrix is constant times the identity: 

1 f°° 

Dij (e) = D e Sij, D e := 2 / £ e [v{t) • u(0)]dt . (2.20) 

[Zir) a J o 

The main result of the paper is the following theorem. 
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Theorem 2.2 Let d > 3 and ifo E T 2 (R rf ) 6e a normalized initial wave function. Let 
if {t) := exp(—itH)ifo solve the Schrodinger equation 11. 1 1 ). Lei 0(x,v) be a Schwarz 
function on M. d x W 1 . For almost all energy e > 0, [|V’o(u)| 2 ](e) finite and let f be the 
solution to the heat equation 

d T f(T, X , e ) = D e A x f(T, X, e) (2.21) 


with the initial condition 


f (0, X, e) :=5(X) \Mv)\ 2 (e) 


for these energies. Then there exist 0 < no(d) < 2 such that for 0 < k < Ko{d) and for 
e and A related by (12.171) . the rescaled Wigner distribution satisfies 

lim J dX j dv 0 (X,v)VW* {x - k . 2t) (X,v) = J dX jdv 0(X,v)f{T,X,e(v)) , (2.22) 

and the limit is uniform on T E [0, To] with any fixed Tq. In d = 3 one can choose 
ko(3) = 1/370. 

Remark 1. The coefficient [|?/>o(u)| 2 ] (e) in the initial condition f(0,X,e) is finite for 
almost all e by using (|2.13l) for h = \ifo\ 2 - 

Remark 2. The total cross section of the collision process (12.181) . 


o'o(e) 


J du a(u,v ) 


e = e(v) , 


(2.23) 


is a function of e = e(v) only. Assuming B( 0) / 0, we see that <ro(e) ~ [1] (e) for small e, 
and cro(e) ~ e^ 1 / 2 for large e. It follows from Lemma 12.II and from standard probability 
arguments that the diffusion constant (12.201) scales as D e ~ e/cro(e) for small e <C 1 
and D e ~ e 2 /cro(e) for large e> 1. If B vanishes at 0 (but (12.51) still holds), then the 
small energy behaviour of ao(e) and D e depends on the rate of vanishing of B at 0 in a 
straighforward way. 

Remark 3. The condition (12.51) is not essential, but the theorem needs to be modified 
if B vanishes on D(0,S), a ball of radius 6 > 0 about the origin. Let 6 > 0 be the 
maximal radius so that D( 0, 5) D supp (B) = 0. In this case the total cross section cro(e) 
is zero for all energy values e < 5 2 / 8, because the diameter of the energy surface E e is 
smaller than the minimal range of B. Therefore the evolution is ballistic for the part of 
the initial wave function that is supported on energy shells e < 6 2 / 8. For the other part 
of the wave function the diffusion equation still holds. 


Fig. □ below shows the three different scales schematically. On the Schrodinger scale 
both time and space are of order 1 in atomic units. On the kinetic scale time and space 
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Diffusive scale: X, T 


Kinetic scale: X, T 


Atomic scale: x, t/ 



--2- K 
Time: A 

Heat equation 




A" 


Boltzmann eq. 


1 

1 

Schrodinger eq. 


Figure 1: Evolution equations on three scales 

are rescaled by A~ 2 . The dynamics is given by the Boltzmann equation characterized 
by finitely many collisions. On the diffusive scale we rescaled the time and space by an 
additional factor X~ K and A _K / 2 respectively. The typical number of collisions is of order 
A 2 t~ A~ K . 

If we assume that the Boltzmann equation holds under all scalings, Theorem 12. 2 1 can 
be easily understood. From the Boltzmann equation (12.161) . the momentum distribution 
develops according to the Markovian generator L e . Therefore, the Boltzmann equation 
(12.161) describes a process that a particle travels with a fixed momentum v up to an 
exponentially distributed random time with average value 0 o(e(r;)) , then it changes 

momentum from v to a new momentum u on the same energy surface X e chosen by the 
probability distribution P(u) = a(u,v)/ao(e(v)). The different energy sectors do not 
interact. Clearly, this process then converges to a Brownian motion in configuration 
space with a diffusion coefficient given by (12.201) and with momentum restricted to a 
fixed energy shell X e . 

Under the assumption that the Boltzmann equation is valid for all time, this argument 
applies in all dimensions. The random Schrodinger evolution, however, is expected 
to be localized for d < 2 even for small coupling constant. Therefore, even though 
the Boltzmann approximation was proved to be valid for d > 2 [16], [10] (it is not 
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valid for d = 1) in the weak coupling limit, it will not be valid for all time in d = 2 . 
It is expected that memory effects and quantum correlations eventually dominate the 
evolution and ruin the Markovian character of the Boltzmann picture. Heuristic ideas 
show this transition happens at an exponentially large time (see, e.g. @71). 

The effects of the quantum correlations and memory are not expected to change the 
Boltzmann picture drastically in d > 3, but one expects corrections to the diffusion 
equation and a transition between different energy shells for n > 2 (see EH]). 

Although Theorem 12.21 is formulated for the Euclidean space, quantum diffusion on 
large scales is not restricted to continuum models. In fact, we prove in |18j an analogous 
result for the Anderson model where the Hamiltonian (|2.1I) is defined on £ 2 (Z rf ), d > 3. 
The kinetic energy is the discrete Laplacian and the random potential is given by 

V u (x) := ^u 7 <5(x — 7 ), x £ Z d , (2.24) 

7 £Z d 

where {u 7 = u 7 (u;) : 7 G Z d } are i.i.d. random variables. 

The discrete model is technically more involved for two reasons. First, the dispersion 
relation of the lattice Laplacian lacks convexity which simplifies several estimates in the 
continuum model. We will address this issue in [El- 

Second, our choice of the potential in the continuum model contains more random¬ 
ness than in the discrete one. Comparing (|2.2[) - (|2.3I) with (12.241) we note that in the 
continuum model both the locations of the obstacles and the coupling constants are ran¬ 
dom variables, while in the Anderson model the locations are deterministic. Formally, 
the potential (12.241) corresponds to (12.21) - ([2.31) with the deterministic choice y 1 = 7 E Z d 
and with B(x) being the lattice delta function 5(x). The random choice of the centers 
in the continuum model simplifies the formula in the calculation of the high moments of 
the random potential. Nevertheless, we will formulate our estimates in a unified setup 
that can be applied directly to the Anderson model in [18] as well. 

2.3 Strategy of the proof 

The above heuristic argument using the Boltzmann equation, besides being misleading 
for d = 2, also masks the difficulties in proving Theorem 12.21 namely that one has to 
follow the full quantum mechanical time evolution through infinitely many collisions. The 
main tool of our proof is to use the Duhamel expansion to decompose the wave function 
into elementary wave functions characterized by their collision histories. We then apply 
two strategies to simplify the expansion: (i) renormalization of the propagator, i.e., 
resumming the two legged subdiagrams; (ii) stopping rules to control recollisions. Apart 
from these two steps, the bulk of our proof is devoted to giving sharp estimates for a 
large class of Feynman graphs. 
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To get an idea, imagine that we expand the solution to the Schrodinger equation 
by using the Duhamel formula repeatedly. This rewrites the solution into a sum of 
elementary wave functions, each of which is characterized by a sequence of collisions with 
random obstacles. When we take the expectation of |v/jj11 2 with respect to a Gaussian 
randomness, we pair the random obstacles by Wick’s theorem and obtain a sum of 
amplitudes of Feynman graphs. In case of a non-Gaussian randomness the higher order 
cumulants are also present due to various recollision patterns (their contribution turns 
out to be negligible, but proving this is rather involved). 

If we take only the Laplacian as the free part in the expansion, even the amplitudes 
of individual graphs diverge in the limit we consider. However, this can be remedied by 
a simple resummation of all two-legged insertions caused by the lowest order self-energy 
contribution (it turns out that higher order corrections to the self-energy do not play 
a role in the scaling limit we consider). The resummation is performed by choosing 
an appropriate reference Hamiltonian Hq for the expansion. After this rearrangement, 
all graphs have a finite amplitude in our scaling limit, and the ladder graphs give the 
leading contribution. 

However, we have to estimate not only individual graphs but the sum of all graph 
amplitudes, which requires beating down the factorial growth of the number of graphs. 
This problem has been addressed in constructive field theory. For field theories with 
bosons, the graphical expansion to infinite order diverges. Borel summability was proven 
by cluster expansion and renormalization group methods [301 E2 EU5]. In fermionic 
theories, the anticommutation relations entail cancellations which lead to analyticity in 
the presence of regulators [HI 123 EH 02]. Our method to control the combinatorial 
growth is completely different: it is by very sharp bounds on the individual graphs. We 
give a classification of arbitrary large graphs, based on counting the number of vertices 
carrying extra oscillatory effects. The number of these vertices is called the degree of 
the graph and it measures the improvement over the standard power counting. For the 
ladder graphs, the degree is zero, for the anti-ladder (i.e., complete crossing) graph it is 
2. For general graphs, the degree is roughly the number of vertices after removing all 
ladder and anti-ladder subgraphs. We thus obtain an extra A c factor (for some c > 0) 
per non-(anti)ladder vertex. This strong improvement is sufficient to beat the growth of 
the combinatorics in the time scale we consider. To our knowledge, nothing like this has 
been done in a graphical expansion before. Improved phase space estimates have been 
used to prove regularity in two-dimensional many-fermion systems, but the improvement 
exponent was fixed independently of the number of vertices [261 [23 [H]. 

For a comparison, the unperturbed Green functions in the perturbation expansion 
for the many—fermion systems and for the random Schrodinger equation are given by 

1 1 
ipo + p 1 — p ’ p 2 — a + irj 
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In the many-fermion case, /j > 0, po E Mp = {-|(2 n + 1 ): n € Z} where (3 ~ T 1 is the 
inverse temperature. In the random Schrodinger case, r] ~ £ _1 . Their T 2 -properties are 
different: 

4 f dp^Po+P 2 “Ml” 2 ~ I log/?|, 

Notice the divergence is more severe for the random Schrodinger equation case. In the 
many-fermion case, there is one /Jo-summation per line of the graph; in the random 
Schrodinger case there are just two overall a-integrals for graphs with arbitrarily many 
lines. 

This paper is organized as follows. In Section [3] we perform the self-energy renormal¬ 
ization, we smooth out the data and restrict the problem to a finite box. The Duhamel 
expansion is introduced in Section |4j In Section [5] we reduce the Main Theorem to 
Theorems 15.11 15.21 and 15.31 The key result is Theorem 15.21 which we prove in the rest 
of this paper. The other two theorems are more technical and they are proven in the 
companion paper m- The Feynman graphs are introduced in Section [7] In Sections 0 
and [9] we reduce all estimates to Theorem 18.41 This theorem is our main technical bound 
on Feynman graphs and it is proven in Section [101 

Since the random potential in our model is given by general i.i.d. random variables, 
the rule for taking the expectation is different from the case of Gaussian random field used 
in [16] . This produces technical difficulties especially for the Anderson model, where, in 
addition to the usual pairing from the Wick theorem, we have to introduce higher order 
partitions of the vertices, called non-trivial lumps. Our continuum model avoids this 
complication due to the additional randomness of the obstacle centers. Nevertheless, we 
present the general proof here so that the key technical results in this paper could be 
applied to the Anderson model as well. For readers interested only in the continuum 
case, we recommend to ignore the non-trivial lumps. 

Universal constants and constants that depend only on the dimension d, on the final 
time To and on the Schwarz norms ||.B||fcfc from (12.61) will be denoted by C and their 
value may vary from line to line. 



3 Preparations 

3.1 Renormalization 

The purpose of this procedure is to include immediate recollisions with the same obstacle 
into the propagator itself. This is also called the renormalization of ” one-particle propa¬ 
gators” or two legged subdiagrams. Without renormalization, these graphs individually 
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are exponentially large (“divergent”), but their sum is finite. Renormalization removes 
this instability and the analysis of the resulting Feynman graphs will become simpler. 
The self-energy operator is given by the multiplication operator in momentum space 


6(p) := 0(e(p)), 0(a) := lim 0 e (a) , 0 £ (a) := 0 e (a,r) (3.1) 

£— > 0 + 


for any r with e(r) = a, where 


0 e (a,r) 


\B(q-r)\ 2 dq 
a — e(q) + is 


(3.2) 


Note that by spherical symmetry of B and e(q), 0 e (a,r) depends only on the length of 
r, therefore 0 e (a) in (|3.1I) is well defined. Clearly 9(p) is spherically symmetric. The 
existence of the limit and related properties of 0 have been proven in [16] using that 
11 B 2 11 2 d, 2 d < oo. Here we summarize the results: 

Lemma 3.1 In d > 3 the following hold: 


|0 e (o:,r) - 0 e (a,r')| < C\ |r| - |r'| | 

(Eq. (3.80) in fl 6 j /) and 

|0 £ (a, r) — 0 £ i(a', r)\ < C(\e — s'\ + \a — a'De -1 / 2 


(3.3) 


(3.4) 


if s > s' > 0 (Eq. (3.68) in )16(). From this latter estimate the existence of the limit 
lim £ ^o+o 0 e (o!, f) follows. Moreover, 0 is Holder continuous 

|0(a) — 0(c/)| < C\a — c /| 1/2 . (3.5) 


Proof. We have only to prove the Holder continuity. For any s and any r, r' with 
a = e(r), a' = e(r') we have 

| 0 (a) — 0 (c/)| < lim |0 e /(a, r) — 0 e (a,r)| + |0 e (a,r) — 0 £ (cc, r')\ 

e'^0+0 

+ | 0 e (o:,r / ) - © e (a',r')| + lim \0 e (a',r') - 0 E :(a',r')\ 

e'^0+0 

< c(e 1 / 2 + | |r| — |/| | + |cc — a'le -1 / 2 ^ . 

By optimizing e and using e(r) = a, we obtain (13.51) . □ 


We have the following estimate on 0(p) and in parallel on 0(e): 
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Lemma 3.2 For any d> 3 there exist universal positive constants ci,c 2 such that 


\m\ < 


C 2 log (p) 


|0(e)| < 


c 2 log(e) 


(3.6) 


(p) ’ " (e) 1 / 2 ’ 

ImQ(e) < —ci min{|e| |e| -1 / 2 } , ImO(p) < — ci min{|p| d ~ 2 , |p| -1 } . (3.7) 

Proof of Lemma, \3.2\ By performing the angular integration, we can write Q e (a,p) 
with e(p) = a as 




(a,p)= [ -—— 5(e) , with 5(e) := [ \B(V2e(4> r ~ 4>))\ 2 d4> , (3-8) 

J 0 a- e + ie J S d -1 

where cj) r is a fixed vector on the unit sphere 5 rf_1 . For small e values 

|5(e)| = 0(l), \VS(e)\=0(e~ 1 / 2 ). 

For large e values, using the regularity of B , 

|S(e)|=0(e-^), |V5(e)|=0(e-i). 

These estimates, in particular, stand behind the proof that lim £ _»o+o 0 £ (a,p) is finite, 
since they guarantee the sufficient decay for large e and the sufficient smoothness around 
the singularity of the denominator in (13.81) . The imaginary part therefore is 

Im#(p) = Im lim^0 £ (a,p) = —7r(2o:)2^ 1 5(a) 

which behaves as ~ ~\p\ d ^ 2 for small p and as ~ — \p\ _1 for large p. The real part of 
Q e (a,r ) is bounded for small a. For large a one splits the integration 


Re 0 e (a, r) = 


r*a+1 


+ 


(2e)z i de 


i-i .7|a-e|>i / a — e + ie 


-5(e) 


After Taylor expanding (2e)z x 5(e) around a» 1, the first term is bounded by 


(2a 


ra +1 

'5(a) / 

Ja-1 


de 


_i a — e + ie 


+ 2 sup 

|e—«|<1 


de 


[(2e) a 1 S(e)] =0(a 1/2 ) , 


and the second term by 


de 


|a—e|>i a e|e 1 / 2 


< 


clog (a) 
(a) 1 / 2 
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If we write 0(e) = H(e) — il(e), where lZ(e) and 1(e) are real functions, and recall 
Im(x + 10) _1 = —7 t5(x), we have 


1(e) = —Im0(e) = 7r J S(e(q) — e)| B(q — r)\ 2 dq 
for any r satisfying e = e(r). □ 


(3.9) 


We rewrite the Hamiltonian as 

H = H 0 + V, 


where 

H 0 :=u(p) :=e(p) + X 2 0(p), V := XV - \ 2 d(p) . (3.10) 

We note that our renormalization is only an approximation to the standard self-consistent 
renormalization given by the solution to the equation 


u>(p) = e(p) + A 2 lim 
£-> 0+0 


\B(p-q)\ 2 dq 
l o(p) - u>(q) + ie 


(3.11) 


Due to our truncation procedure, the definition (|3.1I) is sufficient and is more convenient 
for us. Since e(p) is spherically symmetric, so are 0(p) and co(p). 


The following lemma collects some estimates on the renormalized propagators we 
shall use to prove Theorem 12.21 The proof is fairly simple and will be given in (19]. We 
note that formula (2.8) in [T9| erroneously contain a factor 2 in the exponent of p, the 
correct bound is (13.141) . 


Lemma 3.3 Suppose that X 2 > p > A 2+4k with k < 1/12. Then we have, 


and for 0 < a < 1 


I h(p - q)\dp 
| a — Lo(p) + irj\ 


< 


\h(p-q)\dp 
| a — a ;(p) + ip | 2_ 

I Hp - q)\dp 
| a — e(p) + ip | 


2 —a 


C\\h\\ 2 d,o log A log (a) 

(a) 1/2 {kl - TAF) 

(3.12) 

^c a \\h\\ 2dfi a -2 ! 1- ®) 
-{a)“/ 2 (|«|-y 2H) ’ 

(3.13) 

^ C a \\h\\ 2dfi p-^- a) 

~ (a)»/ 2 (| g | - v/2W> 

(3.14) 


For a = 0 and with h := B 2 , the following more precise estimate holds. There exists a 
constant Co, depending only on finitely many constants Ck from flZBD such that 


f >?\B(p~q)\ 2 dp 
J | a — ZJ(p) — ip\ 2 


< l + C 0 A- 12K [A + |a-a;(g)| 1 / 2 ] . 


(3.15) 
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3.2 Smoothing the initial data and the potential 

In this section we show that it is sufficient to prove the Main Theorem under the as¬ 
sumptions that is a bounded, smooth, compactly supported function and B(p) is 

supported on {|p| < A^ 5 } for any fixed <5 > 0. 

The approximation procedure relies on the following L 2 -continuity property of the 
Wigner transform. If a random wave function is decomposed as if = if± + if 2 , then 


E{0,W$)-E{0,W; 


£ ' 
Ip 1 I 


< c 


( J sup |(9(£, u)|d£^ [||V>1 


+ 


El 


(3.16) 


by a simple Schwarz inequality. (Due to a misprint, the ||V’2 || 2 term was erroneously 
omitted in Section 2.1. of our earlier paper [16].) 

Approximation of the initial data. Let ifo € L 2 and let if n be a sequence of smooth, 
compactly supported functions with H'i/’n — ifo\\ —► 0. We decompose ifo = ifn + idfo — 'f’n)- 
Then 

ff(t) = e~ ltH if n + e~ ltH (ifo - $n) ■ 

Since 

\\e~ ltH {i >0 - $n)W = Wi’O - $n\\ ->• 0 
as n —> oo, uniformly in t, we see that 


lim 


E(0,W^ t) )-E(0,W r 


e \ 


= 0 


uniformly in t (and thus in e), where if n (t) '■= e~ ltH if n is the time evolution of the 
approximated initial data. This means that the approximation procedure is continuous 
on the left hand side of (12.221) . 

Similarly, on the right hand side of (12.221) . we can define / n (T, X , e) to be the solution 
to (12.211) with initial data f n (0,X,e) := 5(X) [IV’nl 2 ] (e). Clearly [IV’nl 2 ] (e) converges to 
[I'i/’ol 2 ] (e) in L 1 (de). Therefore 


f n (T,X,e)^f(T,X,e) (3.17) 

in LfiydX de), uniformly in T. The right hand side of (12.221) is therefore also continuous 
as n —► oo. 

We remark, that if ifo is smooth, e.g. \X p ifo(p)\ < C(p)~ 4d , then a bounded, smooth 
and compactly supported approximant, if n , can be chosen so that [I'i/’nl 2 ] (e) — > [|V’o| 2 ] (e) 
for every e > 0 and then the convergence in (I3.17P also holds in L l (&X) for any e. The 
smoothness of ifo is used only at the point when we explicitly compute the main term 
of the perturbation expansion and identify it with the Boltzmann equation, see m- 
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Propagation estimate. To verify that a truncation is allowed for B, we first need a 
crude propagation estimate. Define the following event for any Z > 0 

Vz :=\co : [ d|/ 4 ,| (y) < Z(k) ,Vfc G ) , 

L J\y-k\<i J 

where |// w | denotes the total variation of the (random) measure n w . For any fixed G 
let iVfc be the number of Poisson points in the ball { x : \x — k\ < 1}. We compute 


E 


'\y-k\<l 


d|/iw| (y) 


d-\-1 


< EN d+1 \v\ d+1 < C d , 


using (12.41) and that is a Poisson random variable with expectation E= vol(unit ball). 
By Markov inequality 


pm < E 


c d 


kez d 


Z d+1 (k) d+1 


= 0(Z 


-d-1\ 


thus we have 


lim P (fi z ) = 1 . 

Z —>oo 

We decompose 

E(d, wj) = E [1 (n z ){6, W?)] + E [1 m(o, wj) 

where !(•) is the characteristic function. On the set Pi c z we use 


(3.18) 


E 


[im(om)}\ < (J S wpm,vm)\ 


^|| 2 P(D |)^0 


(3.19) 


as Z —> oo, uniformly in t (hence in A). 

For lo G Q z we have |V^(x)| < CZ(x) using the decay properties of B. Com¬ 
puting the time derivative of the mean square displacement, we obtain dt{^tiX 2 ^t) = 
[Hi x 2 ]^t) ■ Using [H, x 2 ] = — (V • x + x ■ V) and a Schwarz estimate we have 


dm,x 2 ip t ) < C{^tiX 2 m l2 [ E + \Vuj\At )] 112 


(3.20) 


with E := ('ipt, Hijjt) = (i/jq, Hipo) by energy conservation. We estimate (ip, |V^|^) < 
CZ + CZ(ip,x 2 Tp)H 2 , in particular the energy E is bounded (depending on i/jq and Z). 
From (13.201) we thus have 


(ipt,x 2 ^ t ) < ci(Z, i/’ 0 )t 4 + c 2 {Z,il) Q ) 
on fl z with some constants ci^iZ, ipo). 


(3.21) 
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Approximation of the potential. We define the truncation of B in Fourier space as 
B s (p) := (p))B(p), where ip : M + —» [0,1] is a fixed smooth cutoff function with 

<p(a) = 1 for a < 1/2 and <p(a) = 0 for a > 1. In position space, we have for any M € N, 


\B(x)-B s (x)\ < (x)~ 2d J\(\7 p ) 2d [B(p)[l-p(\ 5 (p))] |dp < C s , M X M (x)- 2d (3.22) 


by using that B is in Schwarz space (12.61) . 

Let 

H 5 := + A [ B s (x - y)d^(y) 

Z JRd 

be the Hamiltonian with the truncated potential. Let 'iff := e~ ltHS if o be the evolution of 
the wave function under the modified Hamiltonian H d . On the set and for t <C A -4 


dtHt-iffW 2 = —2 Im (iff , (H — H 5 )if t ) < C S Z A 13 (^, (z 2 )^) 1/2 < C(Z,5,if 0 )A 5 

by using (13.221) with M = 12. In particular, ift and iff remain close up to time scale 
t ~ A” 2_k , k < 2. This bound, together with the L 2 -continuity of the Wigner transform 
(13.161) guarantees that the truncation of B does not influence the left hand side of (12.221) . 

As for the right hand side of (12.221) . notice that the collision kernel, a(U,V), of the 
momentum jump process (12.181) is restricted to the energy surface e(V) = e(U) = e. 
Therefore U, V are bounded, depending on e, so B(U — V) = B d (U — V) for these 
momenta, if A is sufficiently small. Thus the truncation of B does not influence the right 
hand side of ( 12 . 221 ) . 

Armed with these results, we assume for the rest of the paper that ifo(p) is smooth, 
compactly supported, bounded and B(p) is supported on {|p| < A _<5 } for any fixed 6 > 0. 
We thus extend the convention from the end of Section [2] that general constants denoted 
by C may depend on the truncated version of B and ifo. The same applies to the hidden 
constants in the O(-) and o(-) notations. 


3.3 Restriction to a finite box 

We will reduce the problem to a finite box of size L, L > 1, with periodic boundary 
conditions. In this way, for technical convenience, we avoid the infinite summation in 
(12.31) . Let A l := [— L/2,L/2] d C be a finite torus and let A^ := (Z /L) d be the dual 
lattice. We introduce the notation 

f{x)dp := j- 1 — f(P ) • ( 3 - 23 ) 

peA* 

The integrals and f A , converge to their infinite volume counterparts as L —> oo. Let 
(■, -)l and || ■ ||l denote the scalar product and the norm on L 2 (Al). 
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For any L,M ^$> 1 we consider the random Schrodinger operator 


H' = H’ lm := 


—A + XV' 
2 “ 


M 

K( X ) ■= XX B ( x ~ 


7=1 



y ) d /4 > 


with periodic boundary conditions on A l and := ■ Here { y 7 : 7 = 

are i.i.d. random variables uniformly distributed on Al and {u ; : 7 = 

1,..., M} are i.i.d. variables distributed according to P„ and they are independent of 
the y'. M itself will be random; it is chosen to be an independent Poisson variable with 
expectation |A/,|. The expectation with respect to the joint measure of {M, y 7 ,u 7 } is 
denoted by E'. Sometimes we will use the decomposition 

E ; = E m E® m E® m (3.24) 

referring to the expectation of M, {y 7 } and {u 7 } separately. The parameter L is implicit 
in these notations. In particular, E® M stands for the normalized integral 


1 

|a 7F 



(3.25) 


It is well known that the restriction of the random measure fj, w (see (|2.3p l to the 
box A 1 has the same distribution as In particular, given a realization u of the 
infinite volume random measure we can associate to it the number of points in 
Al (M = M(ui)) and the operator H' u = with random measure We can 

thus realize the random operator H' u on the same probability space as H w . Due to the 
periodic boundary and the nontrivial support of B, the potential of H u and H' u will not 
be the same on A l, but the difference will be negligible far away from the boundary. 

Let xl be a smooth cutoff function, supported on A l, with xl = 1 on A L j 2 and 
|V X l| < CL _1 . Let := XL £~ ltH and let := e~ ltH 'ip' 0 be the two dynamics 

applied to the cutoff initial data 4>' 0 := xlV’o supported on Al- We also define the cutoff 
observable Ol ■= XlC. Clearly 

l i^o E (d L ,w; L{t) ) L = E <a, W* (t) > (3.26) 


for any t. We estimate 

dt\\Mt)-^ L mi <cm-H')Mmi+c\\[H,xL]Mmi • ( 3 . 27 ) 

The second term is bounded by CL~ 1 \\V'i/jL(t)\\ and on Qz it can be estimated by the 
total energy as in (13.201) . With a propagation estimate similar to (13.211) but applied to 
the evolution we easily obtain that the right hand side of (|3.27l) vanishes as 

L —► oo for any t. On the complement set, lo £ Ll c z , we use the uniform bound (13.191) 
and finally let Z —► oo. In summary, we have shown the following 
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Lemma 3.4 Let ip'(t) := e where M is a Poisson random variable with mean 

|A/,|, then 


nmsup iMU, Wfo t) ) 

L—>oo 


whenever f sup„ |C>(£, u)|d£ < oo. □ 


□ 


4 The Duhamel expansion 

We expand the unitary kernel of H = Hq + V (see (|3.10l) ') by the Duhamel formula. Due 
to the restriction to A l, we really work with H' L M = Hq + V', where the renormalized 

free evolution, Hq. is given by u(p) in Fourier space and V' = XV' — A 2 9(p), p G h* L . The 
prime indicates the restriction to A^ and the dependence on L and M. In this section 
we work on A l but we will mostly omit the primes in the notation. 

For any fixed integer A^ > 1 

TV—1 

tpt ■■= e~ ltH lf }0 = ^2 V’nW + ^ N{t ) , (4.1) 

n =0 

with ^ 

j> n (t) :={~i) n [ [dsj]” +1 e- iSn + lHo Ve~ iSnHo V ...Ve- isiHo i/; 0 (4.2) 

Jo 

being the fully expanded terms and 

%(f) := (-*) r dse-V-'^Vi/tN-iis) (4.3) 

Jo 

is the non-fully expanded or error term. We used the shorthand notation 



Since each potential V in (|4.2I) . (|4.3I) is a summation itself, V = —\ 2 9(p) + 

Vy(x) := v 1 B{x — y 7 ), both of these terms in (14.21) and (14.31) are actually big summa¬ 
tions over so-called elementary wave functions, which are characterized by their collision 
history, i.e. by a sequence of obstacles and, occasionally, by 6{p). Denote by T n , n < oo, 
the set of sequences 

7 = (7i)72,---,7n), 7je{l,2,...,M}U{i5} (4.4) 
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and by Wj the associated potential 


Wj := 


AUy 

-A 2 ff(p) 


if 7 G {1,..., M} 

if 7 = i?. 


The tilde refers to the fact that the additional {i?} symbol is also allowed. An element 
76 {1,...,M}U {$} is identified with the potential W=j and it is called potential label if 
7 G {1,..., M}, otherwise it is a &-label. A potential label carries a factor A, a tLlabel 
carries A 2 . 

For any 7 G T n we define the following fully expanded wave function with truncation 

VHy := H )”- 1 Ada,-]? W^e-^W^ ... (4.5) 

Jo 

and without truncation 

H := (-*)" / t [ds i ]" +1 e ~ iSn + lH °W^ n e~ iSnH ° W’ 7 n _ 1 ... e -^H 0 _ ( 46 ) 

Jo 

In the notation the star (*) will always refer to truncated functions. Note that 

V>t ,7 = (~i) f ds e~ l( ' t ~ s)Ho i/j^ s ^ . 

Jo 

Each term (14.61) along the expansion procedure is characterized by its order n and by a 
sequence 7 G T n . We now identify the main terms. 

Denote by Y^ r C T^ the set of non-repetitive sequences that contain only potential 
labels, i.e. 

rr : = {7 = (71, • • • , Ik) ■ Ij € {1, • • ■, M}, 7j r 7 i if H j} ■ 

Let 

« ~ E *.7 

Ter- 

denote the corresponding elementary wave functions. 

The typical number of collisions up to time t is of order X 2 t. To allow us for some 
room, we set 

K := [A - 5 (A 2 f)] , (4.7) 

([ • ] denotes integer part), where 6 = 8(n) > 0 is a small positive number to be fixed 
later on. K will serve as an upper threshold for the number of collisions in the expansion. 
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5 Proof of the Main Theorem 


The proof is divided into three theorems. The first one states that all terms other than 
iffu, 0 < k < K, are negligible. For the precise statement we use the previous notations, 
in particular we recall that the prime indicates the dependence on L , M 


Theorem 5.1 (L 2 -estimate of the error terms) Let t = 0( A 2 K ) and K given by 
62D. If k < Ko(d) and 5 is sufficiently small (depending only on n), then 


lim lim E' 

A^O L —>oo 


$ 


K -1 

E 


k =0 


2 

L 


= 0 . 


In d = 3 dimensions, one can choose kq(3) = o=q. 


The second key theorem gives an explicit formula for the main terms, 7 r . It really 
identihes the so-called ladder diagram as the only contributing term. We introduce the 
notation 


Rri(az, v) 


1 

a — cv(v) + irj 


for the renormalized propagator. 


Theorem 5.2 (Only the ladder diagram contributes) Let n < 33 J^ 36 , e = \ 2 + K / 2 ; 
t = 0(\~ 2 ~ k ), and K given by (14.71) . For a sufficiently small positive 5, for r] = \ 2+K 
and for any 1 < k < K we have 

lim E'II^TIIl = V x (t,k) + o(\^ 6+ ^ k -°^) (5.1) 

lim (<5 l ,E 'Wfjj „)l = W x {t,k,0) + o( \h(- 6 +^ d y-°^) (5.2) 

L—>oo " t,k \ ) 

as A< 1. Here 

\ 9/t) /* /‘ee /» k~\~ 1 

V\(t,k) :^ / J dad(3 e l(a_/3) * J ( d Pi) l^o(pi)| 2 

k+1 k 

x Y[R n (a,p j )R n (^,Pj)Y{\B{p j+1 -p j )\ 2 (5.3) 

3= 1 3=1 

\2k p 2tri p i-oa /• /• k+l ^ __— 

W x (t,k,(D) :=^^jj_jladfl df J (J]d VjjO&Vk+jW^&Vi) 

k-\- 1-r— k 

x J} R n (a,Vj + y )/? r/ (/3, v 3 - y) JJ | B(vj - u j+ i)| 2 . (5.4) 

3=1 3=1 
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We adopt the notation 0(5) in the exponent of A. This always means (const.)5 with 
universal, explicitly computable positive constants that depend on n and that can be 
easily computed from the proof. 

The formula (15.31) is the value of the so-called ladder Feynman graph in the dia¬ 
grammatic expansion of E , ||^ ? r || 2 . We will see in Proposition 17.21 that this expansion 
generates klB k terms, where B k is the number of partitions of a set with k elements (note 
that Bk is almost of order k\). Theorem 15.21 states that only one diagram is relevant; the 
contribution of all the other Feynman graphs is negligible even after summation. The 
extension of (15.11) to the Wigner transform (15.21) is straightforward. Theorem 15.21 is the 
most important step in the proof of the Main Theorem. 

The third theorem identifies the limit of W\(t. k, O) as A —► 0 with the solution 
to the heat equation. We note that the definition (15.41) does not apply literally to the 
free evolution term k = 0; this term is defined separately: 

W x (t, k = 0,0) := J d£du e it£V< e «* 2lm *(«0 £}(£, v)W 0 {ef, v) . (5.5) 


Theorem 5.3 (The ladder diagram converges to the heat equation) Under the 
conditions of Theorem \5.2\ and setting t = A~ 2 ~ K T, we have 


K -1 


lim ^ W\(t, k, O) = j dX j dv 0(X,v)f(T,X,e(v)) , 


(5.6) 


k =0 


where f is the solution to the heat equation (12.211) . 


Proof of the Main Theorem 12.21 using Theorems 15.1115.41 and 15.51 We compute the 
expectation of the rescaled Wigner transform, EWf = EW5 , tested against a Schwarz 
function 

J dX J dv 0(X,v)EWf(X,v) = f dfj du d(£,v)E W?(£,v) = (0,EWf) . 

Combining Lemma 13.41 Theorem 15.11 and the finite box version of the L 2 -continuity 
of the Wigner transform (13.16|1 . it is sufficient to compute the Wigner transform of 
if'(t, K) := Zk=o ^'tk ■ The Wigner transform W^iur) is quadratic in if ', so it contains 
a double sum over k and k' 

w^,,K)= E (■■■)■ 

k,k '=0 

The potential labels are not repeated within ■*/> and if. Moreover, the expectation of a 
single potential in (14.61) is zero. Thus the potential labels in the if and if must pair, in 
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particular taking expectation reduces this double sum to a single sum over k 


K -1 

E'WV«,*) = £ 

k =0 


By using (15.21) and (15.61) together with K = 0(X K 5 ), we obtain Theorem 12.21 □ 

The main result of the present paper is the proof of Theorem 15.21 The proofs of 
Theorem 15.11 and Theorem 15.31 will be given in the companion paper |19j . For the 
reader’s convenience, we summarize below the key ideas of the proof of Theorem 15.11 
from [19]. 

The Duhamel expansion allows for the flexibility that at every new term of the 
expansion we perform the separation into elementary waves, and we can decide 

whether we want to stop (keeping the full propagator as in (Id.31) 1 or we continue to 
expand that term further. This decision will depend on the collision history, 7 . In 
particular, not every error term will be expanded up to the same order N, in some cases 
we may decide to stop the expansion earlier. 

To estimate a non-fully expanded term, we will use the unitarity of the full evolution, 



e d* ^n^-ds 



(5.7) 


Typically we lose a factor of t by using this estimate since the oscillatory character of 
the time integration is lost. We can use this crude estimate only if the fully expanded 
term, ||?/ , *s, 7 l| 2 ) is small, i.e. if 7 represents an atypical collision sequence. Once 7 is 
“sufficiently” atypical, we stop the expansion for that elementary wave function to reduce 
the number and the complexity of the expanded terms. 

There are basically two patterns how a collision history can become atypical; either 
the total number of collisions exceeds the typical number of collisions, 0(A 2 f), or there is 
a recollision. This explains why only the non-repetition terms ?/dT with k < K contribute 
to the main term. 

A recollision is typically penalized by a factor A 2 in the weak coupling environment. 
This is, however, not the case for the immediate repetition of a potential label, 7 j = 
7 j_i_i E {1 The renormalization (|3.10l) compensates for these terms. Up to 

the highest order, the contribution of a sequence with an immediate repetition cancels 
that of the same sequence where the repetition is replaced by a 0-label. Technically, all 
these estimates have to be combined with the key method of the present paper (proof of 
Theorem 15.21) to show that the sum of all klB^ repetition diagrams is sufficiently small 
to compensate for the unitarity estimate (15.71) . 

The result of the current paper (Theorem 15.21) have been fundamentally used in [19] . 
While that paper was already in print, we have improved the possible range of k and 
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the exponent in the error bounds in Theorem 15.21 (compare with Theorem 2.3 [19]) and, 
therefore, several exponents in [19] can be improved. While these improvements are 
minor (and the exponents used in [19] are still correct), we list them in the following for 
the convenience of the readers. Instead of k < 2/(34 d + 39) required in Theorem 2.3 of 
m, the upper bound k < 2/(33d+36) is sufficient. The exponent 1— (^-d+^-)n — 0(5) 
appearing in Theorem 2.3 of m has been improved to ^ — (6 + ^-d)n — 0(5). Using 
the improved exponent in (19.41) instead of | — (-yd + |)k — 0(5) in (5.18) of [T9j, the 
exponent in (4.37) of Proposition 4.6 in [19] can be improved to | — (7 + ^-d^K — 0(5). 
The better estimates on terms (I) and (II) at the end of Sections 4.5 and 4.6 of [19] 
will lead to a somewhat better threshold Ko(d) for k. More precisely, the upper bound 
(4.31) of |49] is changed to k < ( 2 q — 48)/((33d + 36 )q + 108 + 72 d) and (4.32) of [T9] is 
changed to n < (2 q — 72)/((33 d + 36 )q + 108 + 72 d) after correcting a typo —(16 + 2 d)n 
to —(16 + 12 d)n in the exponent of the previous line. In Section 4.6 of [19], the bound 
(4.40) is changed to k < 2/(18g + 33d + 96) and (4.41) is unchanged. These bounds yield 
an explicit k,q (d) > 0 depending on the dimension, so that Theorem 4.1 of [19] holds for 
k < Ko (d) (the explicit upper bound in Theorem 4.1 was a typo, it should have been 
k < Ko(d)). For d = 3, explicitly kq(3) > -^q. 


6 Pairing potential labels 


The wave function 


= H) fc f\- ds j\ k i +1 e. iSk+lH '°V^ k e~ iSkH '° V ^ fe _ 1 ... e- iS2H oV^e- isiH ^' 0 


contains k potential terms with different potential labels. Every term in 


E 


nr ||2 

°t,k II L 




7;7 


has 2k potential terms, and their expectation is 


E' V l V' V’ V', V', V', 

71 72 v lk 7 ( 72 7 k 


( 6 . 1 ) 


Since there is no repetition within 7 and y, and E ; V/ = 0, the expectation in (16.11) is 
nonzero only if there is a complete pairing between 7 and y. Such pairings correspond 
to permutations on /j. = {1,2,...,!’}. We denote by the set of all permutations on 
k elements. 

We recall the AMdentity from Lemma 3.1 of [16] (with a corrected (27T)” 1 factor) 



ie vt r 

2vr 7 k 


da e 


—iat 


k +1 

n 


L = A a - u(pj) + irj 


( 6 . 2 ) 
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for any p > 0. Therefore, we have 


E 


tU.jjt nr ||2 


\2k e 2t v 

~W) 


(7ee fe 


E //., d p d p ^(Pfc+i — Pfc+i) 


(6.3) 


M k 


XE' E II 1 #(Pj+l -Pj)^'fe(i)+1 ~ PaijMiPl^oiPl) 


71>-~,7fc=l j = l 

/k +1 

x / / dad/3 e* (Q -^ J] 


\ “q a “ u(pj) — ir] (3 — uj(pj) +IT] J 


where the summation runs over all ordered fc-tuples (71, • • ■ ,7 k) of {1,2 ,,M} with 
disjoint elements. We compute the expectation, using m 2 = 1 from (12.41) and the factor¬ 
ization of E ; from (13.241) 


k _ 

E n % (Pj+ 1 “ Pj)%3 (p<r(j)+l - Pati)) = p ( a ’ P> P)S(p)B(p) 

3 = 1 


(6.4) 


with 


and 


k 

e(p)— n^^ +i -pj) 


3 = 1 


M k 

2ni yij (p j+ 1 - pj - (p^yj+i - p CT (i))) 

71...->7fc=l i=l 
Ti^Tj 


P(u, p, p) := E M E® M JJ exp 


(6.5) 


( 6 . 6 ) 


We obtain from (16.31) that 


E' 


/II_/,/ nr ||2 


■t,fc ml - x2k dpd P “Pfc +0 

o-eSj. 1 


At 


(6.7) 


xP(cj, p, p)M°(k, p, p)//>o(pi)V>o(pi) 


with 


M°(k, p, p) ' 

(2vr) 2 .7 ,/ R y ^ a - u{pj) - ip f3 - w(pj) + vq ^ 


( 6 . 8 ) 


(with the convention that for j = k + 1 we set the superfluous term B(pj + 1 — pj) := 1). 
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The expectation value in (16.61) can be easily computed to yield a product of delta 
functions since the variables y 7j . are independent. The constraint 7 j 7 ^ 7 j induces only a 
trivial combinatorial factor that becomes irrelevant in the L — > 00 limit. 

If the obstacle centers were deterministic, i.e., y 7j = 7 j, then the constraint 7 j 7 ^ 7 ^ 
has more serious consequences. This is the case for the lattice Anderson model, where 
the summation in (16.61) extends to all 7 j E Tfi with 7 j 7 ^ 77 and the momentum variables 
are on the dual torus, T rf := [—|, ^]- Due to the constraint 7 * 7 j-, the formula (16.61) is 
not a simple product of delta functions and we have to use a connected graph expansion 
that is well-known from statistical physics. 

Let A n be the set of partitions of Ik := {1,2,... , A:}, i.e. A = { A M : /r E -1(A)} E Ak 
if U = h; and the elements of A are disjoint and non-empty. The sets in the 
partition are labelled by the index set /(A) and let m{ A) = |I(A)| denote the number 
of elements in A. The elements of the partition A will be called lumps. A lump is trivial 
if it has only one element. The trivial partition, where every lump is trivial, is denoted 
by A 0 . 


Lemma 6.1 (i) [Continuum model] For any fixed L, k and M, k < M, and any fixed 
momenta qj E A* L , 


T7\ (g)M 

% 


M k 

e n>p 

71,..., 7fc =l j= 1 


^iQjV'yj 


Ml 

Al|^(M - fc)! 


e n “(I* 4 - 

AgA/c ^g/(A) 


E ( 6 - 9 ) 


with c(l) = 1, c(n) = 0 for any n > 2. 

(ii) [Lattice model] For any fixed k and momenta qj E T c/ , 


E n«p 

71.■■■>7 k^ zd 1 
7i7^7j 


2'xiqj^fj 


e n c (i^i)^( e «<) 

AGA& i/G/(A) -^GAi/ 


( 6 . 10 ) 


with 

c(n) := ^ (-l) |r| 

rcK„ 
r connected 

where K n denotes the complete graph on n vertices and |T| denotes the number of edges 
in the subgraph T. The following estimate holds for n > 2 


c(n)| < n n 2 . 


( 6 . 11 ) 


Remark 1. Recall that M is a Poisson random variable with expectation |A/-J. There¬ 
fore, apart from the prefactor that converges to 1 almost surely as L —» 00 , the right 
hand side of B is simply Y\ k j=1 S(qj). With an obvious choice of c(n) we write it in 
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the same form as (16.101) . In this way the continuum and lattice models can be treated 
simultaneously. The explicit form of c(n) will not be needed. The arguments in the 
sequel will use only the bound (16.111) that is valid for both choices of c(n). 

Remark 2. Analogous formulas hold if the natural index set Ik = {1,2, is 

replaced by an arbitrary finite set S. In this case, the summation on the right hand side 
of (16.91) - (16.101) is over all partitions of S. The set of these partitions is denoted by A(S). 


Proof of Lemma \6. 11 Part (i) is straightforward from the definition of E y . For part 
ii) we use the connected graph expansion 


k 


JJ (1 ^7i ,7? ) 

i+j =1 


e n ^)’ 

AG»4.fc k'EZ(A) 


where 

6 c (A) = c(\A\) 5^, 

e,£'eA 

is the Ursell coefficients of the hard-core lattice gas (see: e.g. m)- Therefore 

L.h.s of dnunj) = E E II n 

7 i,..., 7 fc 6 Z d AeAk vei(A) V £,£'&A V 


= E II [cd^iwE®) 

AGA/j j/£/(A) 


□ 


We will use the identity (16.91) to express P in (16.61) as a linear combination of products 
of delta functions of the momenta and insert it into (|6.7I) . After the limit L —> oo, each 
term in the summation will be expressed by a Feynman graph. The precise 

definitions will be given in the next Section. 


7 Graphical representation 

Traditionally, the Feynman graphs consist of interaction vertices and particle lines among 
them. In case of Gaussian random potentials, the interaction vertices are paired accord¬ 
ing to the Wick theorem m ■ For non-Gaussian randomness, the non-vanishing higher 
order cumulants correspond to joining several vertices [21]. In our case, the appearance 
of the non-trivial subsets are due to selecting the non-repetition sequences. This re¬ 
quires us to define Feynman graphs in a more general setup than usual. In this section 
we introduce the necessary graphical representation in full generality and we will define 
the value of a Feynman graph, V°(A, a), with permutation a and partition A in (17.211) . 
The final result of this section is given in Proposition 17.21 at the end. 
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7.1 Circle graphs and their values 

We start with an oriented circle graph with two distinguished vertices, denoted by 0, 
0*. The number of vertices is N. The vertex set is V, the set of oriented edges is £(V). 
For v G V we use the notation v — 1 and v + 1 for the vertex right before and after v 
in the circular ordering. We also denote e v - = [y — l,u) and e„+ = (v,v + 1) the edge 
right before and after the vertex v, respectively. In particular ep, +1 )_ = e v +. For each 
e G £(V), we introduce a momentum w e and a real number a e associated to this edge. 
The collection of all momenta is denoted by w = {w e : e G £(V)} and dw = <g) e d w e is 
the Lebesgue measure. We sometimes use the notation v ~ e to indicate that an edge e 
is adjacent to a vertex v. 

Let P = {Pfj, : p G 7} be a partition of the set V \ {0, 0*} 

V\{0,0*}= U 

(all P/, nonempty and pairwise disjoint) where I = 7(P) is the index set to label the sets 
in the partition. Let m(P) := |7(P)|. The sets are called P -lumps or just lumps. If 
two elements v, v' G V \ {0, 0*} belong to the same lump within a partition P, we denote 
it by v = v' (mod P). We assign a variable, u M G M d , p G 7(P), to each lump. We call 
them auxiliary momenta ; they will be needed for a technical reason. We always assume 
that the auxiliary momenta add up to 0 

E u » = 0 • ( 7J ) 

The set of all partitions of the vertex set V\{0,0*} is denoted by Vy. For any P C V, 
we let 

L+{P) := {(u,u + l) G Oy) : v + lgP, veP} 

denote the set of edges that go out of P, with respect to the orientation of the circle 
graph, and similarly L_(P) denote the set of edges that go into P. We set L(P) := 
L+{P) U L-{P). 

For any £ G we define the following product of delta functions 

A(P,w,u) :=<$(£+ E ±We ) n < E ±w e - u^j , (7.2) 

eGl/±({0*}) ^E/(P) 

where u := {u^ : p G 7(P)} G W l is a set of auxiliary momenta. The sign ± indicates 
that momenta w e is added or subtracted depending whether the edge e is outgoing or 
incoming, respectively. The function A(- • •) = Ag(- • •) depends on £, but we will mostly 
omit this fact from the notation. All estimates will be uniform in £. 
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Summing up all arguments of these delta functions and using (17.11) we see that these 
delta functions force the two momenta corresponding to the two edges adjacent to 0 to 
differ by £: w e — w e > = £ for e € L + ({ 0 }), e' € L_({ 0 }). 

As a motivation for these definitions, we mention that the lumps naturally arise from 
the connected graph formula (Lemma 16.1|) . According to this formula, the Kirchoff Law 
must be satisfied for all lumps, i.e. the incoming and outgoing momenta must sum up 
to zero. This fact would be described by the delta functions (17.21) with all u^ = 0. In 
certain recollision terms, however, the non-repetition condition leading to Lemma 16.II is 
not fully satisfied and the Kirchoff Law breaks down for a few lumps. The nontrivial 
auxiliary momenta will bookkeep this deviation from the Kirchoff Law (see fl9] for more 
details). Finally, the shift by £ at the vertex 0* in (17.21) will be used when computing 
the Wigner transform in Fourier representation (12.101) . 

For each subset Q C V \ {0,0*}, we define 

w) : = II l^oKOI JJ I B{w ev _ - w £v+ )| - w ev+ y 2d . (7.3) 

e~0 uGV\{O,O*}\0 

In our application, the subset Q collects those vertices, where the original potential decay 
I B(w ev _ — w ev+ ) | could not be explicitly kept along the estimates and this will happen 
only at a few places; the size of Q will be at most 8. For the purpose of this paper, 
i.e. for the proof of Theorem 15.21 we will need only Q = 0, but for the analysis of the 
repetition terms in m we need the more general definition. 

Due to the support properties of B and 'i/jo, we will see that all intermediate momenta 
w e satisfy \w e \ < N\~ s . The maximal number of vertices in our graphs will be N < 
2A'+2 = 0( A -K ~ 5 ), therefore all intermediate momenta will be smaller than £ := A~ K ~ 3S . 
This justifies to define the restricted Lebesgue measures 

d/i(w ) := l(|tc| < C)du> , £ := A~ K ~ 35 , d^(w) := <S> e d/j,(w e ) . (7-4) 


Moreover, each auxiliary momenta, will always be a sum (or difference) of different 
w e momenta (see (E2D), therefore each of them always satisfies \u^\ < 0(A 2k 45 ). 
We will often take the supremum of all possible auxiliary momenta and sup u is always 
considered subject to this bound. 

With these notations, we define, for any P e Vy and ^ = 0,1,2,..., the A-value of 
the partition 


E g ( P,u,a) := A ^ -2 


sup 

g : \G\<g - 


d^(w 


n 

ee£(V) 


a e - Lo(w e ) + ir}\ 


A(P, w, u)A/g(w) 


(7.5) 

The prefactor A N ~ 2 is due to the fact that in the applications all but the two distinguished 
vertices, {0, 0*}, will carry a factor A. The A-value depends also on the parameters A, r/, 
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but we will not specify them in the notation. In the applications, the regularization r/ 
will be mostly chosen as rj = \ 2+K . 

We will also need a slight modification of these definitions, indicated by a lower star 
in the notation: 


E* g { P,u, cc) := A 


N—2 


sup 

g : \g\<g- 


d M w ) n 


ee£(V) 

e^L({0*}) 


\a e - u(w e ) +ir]\ 


A(P, w, u)A/g(w) 


(7-8) 

The only difference is that the denominators carrying the momenta associated to edges 
that are adjacent to 0* are not present in E* g . We call E* g the truncation of E g . We 
will see that Feynman diagrams arising from the perturbation expansion can naturally 
be estimated by quantities of the form (17.51) or (17.61) . 


7.2 Feynman graphs 

We apply this general setup to the following situation that we will call Feynman graph. 
Every quantity in our perturbation expansion will be expressed by values of Feynman 
graphs that are defined below. 

For experts we mention that our Feynman graphs differ from those that one typ¬ 
ically obtains after averaging over a Gaussian disorder. In the latter, potential lines 
never appear as external lines but only as pairing lines, and one can identify vertices 
connected by pairing lines so that the graph becomes four -valent. In our case, the graph 
is still trivalent and has external potential lines, with a corresponding dependence on 
momentum variables u. Also, averaging the disorder will not simply pair up lines but 
can also join more than two potential lines, which correspond to the higher moments. 

Consider the cyclically ordered set V n ,n' ■= {0, 1 , 2, ... , n, 0 *,n' , n' — 1 ,..., 1 } and 
view this as the vertex set of an oriented circle graph on N = n + n' + 2 vertices. We 
set I n := {1, 2 ,... n} and I n ’ := { 1 , 2 ,... , n'} and the vertex set can be identified with 
V n ,n> =/„u7 n /U{0,0*}. 

The set of edges jC(V n y) is partitioned into C(V ntn /) = C U C such that C contains 
the edges between I n U {0,0*} and C contains the edges between I n > U {0,0*}. 

Let V n . n ' be the set of all partitions P on the set I n U I n >. The lumps of a partition 
containing only one vertex will be called single lumps. The vertices 0 and 0* are not 
part of the partitions hence they will not be considered single lumps. Let G = G(P) 
be the set of edges that go into a single lump and let g(P) := |G(P)| be its cardinality. 
In case of n = n', we will use the shorter notation V n = V njn , V n = V n . n etc. The 
Feynman graphs arising from the non-repetition terms will always have n = n' and no 
single lumps, g{ P) = 0, but the more general definition will be needed for the repetition 
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terms in |19] . We remark that even in ESI we will always have 


n — n'\ < <?(P) < 4, to, n! < K . 


(7.7) 


We also introduce a function Q that will represent the momentum dependence of 
the observable. In our estimates, we will always bound Q in supremum norm; no decay 
or smoothness will be necessary. We will need extra conditions on the observable only 
to evaluate the ladder in the proof of Theorem 15.31 (see m for details). Since Q will 
always appear linearly in our formulae, we can assume, for convenience, that ||Q||oo < 1- 
General Q can be accommodated by a multiplicative factor ||Q||oo in the final estimate 
but it will not be carried along the proofs. 

We define the following function to collect all potential terms: 


n [-»(”«)] n 

-0{w e )\ 

][ B(w e -w e+ 1 ) B(w e -w e + 1 ) 

ee£nG ee£nG 


e£C\G 

e/0* 

eeC\G 

e'/'O 


tol l-j 

8 s 

CT) 

o 

1 


(7.8) 

X ^(w eo+ )^(w eo _)Q 

+ »e 0 , + ) 



with w := {w e : e £ £ U £} and recalling that for any e £ £(V), the edge e + 1 denotes 
the edge succeeding e in the circular ordering. 

The delta function A(P,w, u = 0) ensures that the two momenta adjacent to each 
single lump coincide. This holds even for £ ^ 0, recall that 0, 0* are not considered 
lumps. Therefore the distribution _Ad(w)A(P, w, u = 0) is supported on the regime 
with \w e \ < ( for all momenta w e , thanks to the support properties of i/jq, B and to the 
control on the number of terms, to, 7 ?/ < K (see Section [7.II) . In particular 

Af (w)A(P, w, u = 0)dw = A4(w)A(P, w, u = 0)d/x(w) . (7.9) 

Using the boundedness of 6 and |'0o(^)| A C(w)~ 10d , we easily obtain 

|A4(w)|A(P, w, u = 0)dw < C 9(p) A/g(w)A(P, w, u = 0)d/r(w) , (7-10) 


where Q is the set of single lumps and <?(P) = \Q\, since the delta function also guarantees 
that there is no additional decay at the vertices v £ Q in Afg(w) (the last product in 
m is a constant). 

Let cc, (3 € R, P £ V n y and 


V(P,a,0) := A"+"' +! ' <p) f dwj| 




a - uj(w e ) - it] aa _/3 _ u;(w e ) + irj 
e&L ee£ 


(7.11) 


xA(P,w,u = 0)AI(w) . 
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Thanks to (17.91) . the integration measure could be changed to dyu(w). The truncated 
version, P*(P, a, /3), is defined analogously but the a and (3 denominators that correspond 
to e G L({0*}) are removed. 

We set Y := A -100 and define 

p 2trj r rY 

V W (P) := j^y J J dad(3 e^-^V w (P,a,/3) (7.12) 

and 

QitTJ P PY 

E(*)g(P,u) ■= J J dadp Etf g (P,u,a) , (7.13) 

where a in E^ g (P, u, a) is defined as a e = a for e G C and a e := j3 for e G £. The 
notation (*) indicates the same formulas with and without truncation. We will call these 
numbers the V-value and E-value of the partition P, or sometimes, of the corresponding 
Feynman graph. Strictly speaking, the V- and the E- values depend on £ through 
A = A^. When this dependence is important, we will make it explicit in the notation, 
e.g. V = V£. The V -value depends on the choice of Q as well. When necessary, the 
notations V^(P;Q) will indicate this fact. 

Clearly, by using (17.101) . 

|V W (P)| < (CA)s( p ) £7 Wfl (P,u = 0) (7.14) 

with g = <?(P). We will use the notation Eu^ g ( P) := E^ g ( P, u = 0). 

As we will see in (16.21) . for the graphical representation of the Duhamel expansion 
we will really need 

p 2trj t■ /■ 

V£)(P) := J y K dad/le it (“-^P w (P,«,/l) , (7.15) 

i.e. a version of 14* j(P) with unrestricted dad/3 integrations. (The circle superscript 
in V° will refer to the unrestricted version of V). However, the difference between the 
restricted and unrestricted 14-values are negligible even when we sum them up for all 
partitions: 


Lemma 7.1 Assuming that g > A 2+4k and (17.71) . we have 

E k*)( p )-VS)( P ) = 0(A 5(nW) ) . 


(7.16) 


The same result holds if (P) were defined by restricting the a, /?-integral to any 
domain that contains [— Y, Y] x [— Y, Y - ]. 
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Proof. First we consider the case n, n' > 1. To estimate the difference, we consider 
the integration domain where either |a| > Y or |/3| > Y. We assume, for definiteness, 
that |a| > Y. and we estimate all a denominators trivially, 

1 < C_ 

\a — uJ(w e ) — irj\ ~ (a) ’ 

by using that \uo{w e )\ < + 0(A 2 ) < + 0(A 2 ) on the support of dfi(w e ). Then we 

estimate all but the last /3-denominators in (17.111) trivially by ri^ 1 (P/Y)~ 1 . Thus all w e 
integrations are trivial except the last one where we use (13.12[) . Thanks to the bounds 
fy(w)\, \B(w)\ < C(w)- 10d , one easily obtains that 


\V(P,a,0)\ < 


( CA )n + n' + g( p )|iogA| jpg^) 

(a) n+1 r] n ' ((3/Y) n ' (/3) 1 / 2 


Therefore, we have 


[ dp( 

JM. U 


n«i>n 


da \V(P,a,(3)\ < 


(C A) 


n+n'+g{ P) | 


v n' Y 


n-i-25 


= 0(A 6( n+n*) ^ 


by using (H2D- Similar bound holds for the truncated values, V*. Thus 


y w (p)-'F ( : ) (p) 


0(A 6 ( n+n ')) . 


(7.17) 


Since the total number of partitions, \V n)n '\, is smaller than (n + n') n+n ' and in our 
applications n, n' < K <C A _K ^ 25 , we see that the restriction of the a,/3-integral to any 
domain that contains [— Y, Y] x [—Y,Y\ has a negligible effect of order 0(X 5 ^ n+n '^) even 
after summing up all V-values. 

If the condition n,n' > 1 is not satisfied, say the number of a-denominators is one 
(n = 0), then we will not introduce the auxiliary variable a as in (16.21) . because the f da 
integral would be logarithmically divergent after taking the absolute value. In this case, 
we use the definition 

V°(P) :=p f d/3 e~ it/3 V (P, /3) 

2vr Jr 


with 


p°(p,/3) : = A n+n,+9(p) j (n 


d w, 


eecuc 


n 


/3 - Lu(w e ) + if] 


Xe i M«')A(P, w , u = 0)M(w) , 


directly, instead of V° given in (17.151) . Similar modifications hold for the other cases 
(' n' = 0, n > 1 and n' = n = 0) as well. In particular, in our expansion (including all 
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Figure 2: Vertex and edge labels 

the cases in m) only a few such graphs may appear due to \n — n'\ < 2. The estimates 
leading to (17.161) in these cases are similar but much easier than in the general case and 
they are left to the reader (the same estimates were covered in [16] as well, without the 
renormalization of the dispersion relation). □ 

Sometimes we will use a numerical labelling of the edges (see Fig. [2]). In this case, 
we label the edge between (j — 1 ,j) by ej. the edge between (j , j — 1) by e--. At the 
special vertices 0,0* we denote the edges as follows: e n+ \ := (n,0*), := (O*,^), 

ei = (0,1) and er := (1,0). Therefore the edge set C = C(V n y) is identified with the 
index set I n+ \ U I n '+i and we set 

Pj := w ej , pj := w e -. (7.18) 

These two notations will sometimes be used in parallel: the p, p notation is preferred 
when distinction is needed between momenta on C and C edges and the w notation is 
used when no such distinction is necessary. Note that we always have 

Pi~Pi = €- (7.19) 

7.3 Non-repetition Feynman graphs 

A partition P G V n of I n U7 n is called even if for any £ P we have |P M n/ n | = |P M nI n |. 
In particular, in an even partition there are no single lumps, G(P) = 0. 

Let & n be the set of permutations on I n and let id be the identity permutation. Note 
that A G An and a G & n , uniquely determine an even partition in P(A,cr) G V n . by 
/(P) := 1(A) and := U <t(A m ). 

Conversely, given an even partition P G V n , we can define its projection onto I n , 
A := 7r(P) G A n , by 1(A) := I( P) and A M := P ^ H I n . We let 

©n(P) :={<tG6„ : P(7t(P),<t) = P} 
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be the set of permutations that are compatible with a given even partition P. In other 
words, a £ 6 n (P) if for each i £ I n the pair (i,a(i)) belongs to the same P-lump. 
Clearly 

l©n(P)|= n (^) ! = n 1^1 !. (7.20) 

A*G/(P) /xe/(ir(A)) 

We will use the notation 

V^(A,a,Q) := V W (P(A,a); Q) (7.21) 

and similarly for E(*) g and . In the proofs, Q will be omitted. We also introduce 

C (A) ]1 C (IAI). < 7 ' 22 ) 

u£l(A) 

With these notations we can state the representation of the non-repetition terms as a 
summation over Feynman diagrams: 


Proposition 7.2 With Q = 1 and £ = 0 we have 

lim E'll^ni = E E c(A)P e ° =0 (A,a,Q = l) (7.23) 

L— KX) ’ z —' z ^ 

o-gSfc A eAk 

and with Q%(v) := 0(£,v) we have 

'(d L ,Wfo) L = E E C ( A ) / de^(A,a,Q ? ) . (7.24) 

<rE(5fc A^Ak 


L 


Proof of Proposition | 7. 4 We insert (16.611 and (16.911 into (16.711 and we take the limit 
-> oo. We use that 

Ml 


E 


M 


1 


l\A L \ k (M-k)ll 

for any fixed k. We also replace every Riemann sum (j3.23f) with integrals and we use 
V’o —> V’o- By recalling (17.191) and by choosing Q = 1 in the definition (17.81) . we obtain 
(17.231) . The proof of (17.241) is identical. □. 


Proof of Theorem 1 5. 71 We will prove only (15.11) : the proof of (15.21) is analogous. 
Starting with (17.23p . we notice that the graph with the trivial partition Ao and with the 
identity permutation on I gives the main term in Theorem 15.21 since 


V\(t,k) = V^ =0 (A 0 ,id) . 


This graph is called the ladder graph (Fig. [3]). All other graphs will be negligible. 
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Ladder graph 


Figure 3: Ladder graph 

We first replace V°(- ■ ■) with V(- ■ ■); the error is negligible by Lemma f7.ll In Section 

[8] we then estimate V(A.,a) for the trivial partition A = Ao, where every lump has one 
element. The result is Proposition 18.61 In this case we set V{a) := V(Ao, cr). In Section 

[9] we treat the general case A ^ Ao- The final result of this section is Proposition 19.21 

The proof of both propositions rely on Theorem 18.41 that is the core of our method. Its 
proof is given separately in Section [101 Finally, the proof of Theorem 15.21 follows from 
Proposition 19.21 together with (17.231) . (17.161) and (|6.11D . □ 

We remark that the E- and F-values of the partitions depend on the parameters 
A, t, £, £, k and g\ a fact that is not explicitly included in the notation. In Sections 0 [9] 
and m we will always assume the following relations 

7] = A 2+K , t = X~ 2 ~ K T, T E [0, T 0 ], K = [A -5 (A 2 f)], k < K, ( = X~ K ~ 35 , g< 8 

(7-25) 

with a sufficiently small positive 5 > 0 that is independent of A but depends on n. All 
estimates will be uniform in £ and in T £ [0,To]. We mention that for the proof of 
Theorem 15.21 we need only g = 0, but the more general case is used in [19]. 


8 Estimates on Feynman graphs without non¬ 
trivial lumps 


We use the letters pj. pj, j E Ik+i f° r the momenta variables (see the convention at 
the end of Section [7.21) and /(Ao) = Ik for the index set of the trivial partition. In the 
following sections we always assume Q = 1. 

We introduce the restricted version of M° (see (I6.8f) ]l as 


M(k, p,p) := 


e 2tr> 


Y 


dad/3 


-Y 


k -\-1 

n 


3 =1 


B(Pj+i ~Pj) B(Pj +1 ~Pj) 
a — ZJ(pj) — ip (3 — Lo(pj) + ir] 


( 8 . 1 ) 


and we also define the trivial estimate of M as 


N(k, p,p) := 


e 2tri 


rY fk +1 

l^dadpl n 


\B{p j+ i-Pj)\ | B(p j+ 1 -pj )| 


\j =i 


\a - u(pj) - ir]\ \P - u(pj) + irj\ 


( 8 . 2 ) 

















The truncated versions of these quantities, denoted by M*(A;,p,p) and N*(k, p, p), are 
defined by removing the (.k + l)-th a and (5 denominators from the definitions (18.11) . 
(18.21) but keeping all numerators and all other denominators. 

From the definitions (I7.12|) . (17.211) and V(a) = V(Ao,<t) = V(Ao, cr,Q = 1), we 
obtain 

V(*)(v) = A 2A II dpdpM(*)(A;,p,p)A f (<T,p,p,u = 0)^o(pi)^o(pi) , (8.3) 

E(*)(cr,u) = \ 2k JJ dpdp A w (A;,p,p)A € (cr,p,p,u)V! o (p 1 )'0o(Pi) (8.4) 

with 

k 

Ag(cr,p,p,u) := 5(p k+ 1 - Pk+l + - Pe - (Pa(e )+1 ~Pa(i)) ~ U t ) • 

£=1 

Clearly | V(*)(cr)| < E^ g=0 (a, u = 0) for any £ (see (1 7.14 [) ). 

We introduce a convenient notation. For any (k + 1) x (k + 1) matrix M and for any 
vector of momenta p = (pi,.. .p k + i), we let M p denote the following (A; + l)-vector of 
momenta 

k -\-1 fc+1 

M P := ( M 2• • •) • ( 8 - 5 ) 

i=i 1=1 

Furthermore, we introduce the vector v = (vi,..., v k +\ ) as 

vn := £ + + U2 + ... + ui- 1, for all K = 1, 2,... , k + 1 . (8-6) 

Note that v k +i = £ by (17. 1 [ ) . 

Given a permutation cr £ ©/,, we define a (A; + 1) x {k + 1) matrix M = M(a ) as 
follows 

T 1 if a(j — 1) < i < a(j) 

Mij{a) := < -1 if a(j) < i < a(j - 1) (8.7) 

[ 0 otherwise, 

where, by definition, a is the extension of cr to a permutation of {0,1,..., k + 1} by 
cr(0) := 0 and a(k + 1) := k + 1. In particular [M p]i = pi, [Mp] k+ i = p k +i- 
It is easy to check that 

k+l 

A 5 (ct, p, p, u) = Pj - [Mp] j + [Mvjjj , (8.8) 

3 = 1 

in other words, the matrix M encodes the dependence of the p-momenta on the p- 
momenta and the u-momenta. This rule is transparent in the graphical representation 
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Figure 4: Domain of dependencies of the momenta 

of the Feynman graph: the momentum pj appears in those pf s which fall into its ” domain 
of dependence”, i.e. the section between the image of the two endpoints of pj , and the 
sign depends on the ordering of these images (Fig. 01). 

Definition 8.1 A matrix M with entries 0,+l or —1 is called tower matrix if in 
each column the non-zero entries are consecutive and identical. The collection of these 
consecutive 1 or — 1 entries are called the tower of that column. 

By construction, the matrix M ( a ) is a tower matrix. 

Proposition 8.2 For any permutation a £ &k the matrix M(a) is 

(i) invertible; 

(ii) totally unimodular, i.e. each subdeterminant is 0 or ±1. 

Proof. The invertibility follows from the fact that p and p play a symmetric role in 
(|8.8|) if v = 0, £ = 0, in particular M(ct) -1 = M(a ~ 1 ). It is easy to prove by induction 
on the size of the matrix that any tower matrix is totally unimodular. □. 

The following definition is crucial. It establishes the necessary concepts to measure 
the complexity of a permutation. 

Definition 8.3 (Valley, peak, slope and ladder) Given a permutation a £ (3/ c let a 
be its extension. A point (j, <r{j)), j £ Ik = { 1,2,..., k}, on the graph of a is called peak 
if a(j) < min{cr(j — 1 ),cr(j + 1)}, it is called valley if a(j) > max{cf(j — l),d (j + 1)}. 
Furthermore, if a(j) — 1 6 {a(j — 1 ),&{j + 1)} and (j,cr(j)) is not a valley, then the 
point ( j,cr(j)), j £ Ik, is called ladder. Finally, a point (j, cr(j)), j £ Ik, on the graph 
of a is called slope if it is not a peak, valley or ladder. 

Let I = {1, 2,..., k + 1} denote the set of row indices of M. This set is partitioned 
into five disjoint subsets, I = I p U I v U Ig U I s U Ii as t, such that Ii ast := {k + 1} is the last 
index, and i £ I p ,I Vl Ii or I s depending on whether (c r“ 1 (z),*) is a peak, valley, ladder 
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or slope, respectively. The cardinalities of these sets are denoted by p := \I p \, v := \I V \, 
t := \Ig\ and s := |J S |, and, if necessary, we indicate the dependence on a as p = p{cr), 
etc. We define the degree of the permutation a as 

deg(cr) := k — £(a) . (8-9) 

A maximal collection of consecutive ladder indices, i + 1,..., i + b e Ig, is called a 
ladder of length b. The index i is called the top index of a ladder. The bottom 
index of a ladder is defined to be i + b or i + b+l, depending on whether |d -1 (i + b + 
1) — a~ l (i + 6)| / 1 or \d~ l (i + b + 1) — cr _1 (z + b)\ = 1, respectively. The set of bottom 
and top indices are denoted by R, It and It C {0,1,... , k — 1}, R C {1, 2,... , k + 1}. 
Note that the top index of a ladder never belongs to Ig and it may be 0. The bottom 
index is either a ladder, a valley or k + 1. 

Remarks: (i) The terminology of peak, valley, slope, ladder comes from the graph of 
the permutation d drawn in a coordinate system where the axis of the dependent variable, 
cr(j), is oriented downward (see Fig. [5]). It immediately follows from the definition of 
the extension d that the number of peaks and valleys are the same, p(cr) = v(a). 

(ii) The nonzero entries in the matrix M(a) follow the same geometric pattern as 
the graph: each downward segment of the graph corresponds to a column with a few 
consecutive l’s, upward segments correspond to columns with (—l)’s. On Fig. [5]we also 
pictured the towers of M (cr) drawn inside the graph of a. 

(iii) Because our choice of orientation of the vertical axis follows the convention of 
labelling rows of a matrix, a peak is a local minimum of j —> a(j). We fix the convention 
that the notions “higher” or “lower” for objects related to the vertical axis (e.g. row 
indices) always refer to the graphical picture. In particular the “bottom” or the “lowest 
element” of a tower is located in the row with the highest index. 

Also, a point on the graph of the function j —> a(j) is traditionally denoted by 
(j, cr(j)), where the first coordinate j runs on the horizontal axis, while in the labelling 
of the (z,j)-matrix element Mij of a matrix M the first coordinate i labels rows, i.e. it 
runs vertically. To avoid confusion, we will always specify whether a double index (i, j) 
refers to a point on the graph of cr or a matrix element. 

(iv) We note that for the special case of the identity permutation a = id we have 
I p = I s = I v = 0, and Ig = { 1,2,... , k}. In particular, deg(id) = 0 and deg(cr) > 2 for 
any other permutation cr 7 ^ id. 

An example is shown on Fig. [5] with k = 8 . The matrix corresponding to the permu- 
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Figure 5: Graph of a permutation with the towers 
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tation on this figure is the following (zero entries are left empty) 




1 


1 

1 

1 


-1 1 
-1 


\ 1 l 

2 e 

3 p 

4 t 


M(a) : 


1 5 s 

1 6 t 

1 7 v 

1 8 s 


( 8 . 10 ) 


1 -1 
1 -1 


V 


1/ 9 (last) 


The numbers on the right indicate the column indices and the letters show whether 
it is peak/slope/valley/ladder or last. In this case I p = {3}, I v = {7}, I s = {5,8}. 
h = (1) 2,4, 6 }, Ii as t = {9}, It = {0,3, 5 }, I b = { 2 ,4, 7} and deg(u) = 4. There are three 
ladders, two of them of length one and one is of length two. 

Now we are ready to estimate | V’(<t)| < E g= o(a, u = 0). The following theorem shows 
that the degree of the permutation, deg(cr), measures the size of V(a). The proof is the 
key step in our method and it is given in Section [TOj 

Theorem 8.4 Assume (17.251) with k < and let a E ©&. Then the E-value of the 
graph of the trivial partition with permutation a is estimated by 



( 8 . 11 ) 


if A< 1 . 


This theorem is complemented by the following combinatorial lemma. 

Lemma 8.5 Let k < K = 0(X~ K ~ 5 ), D > 0 integer, and let 7 > k + <5 be fixed. Then 



( 8 . 12 ) 


cG6 

deg (a)>D 


for A< 1 . 


Since deg(u) > 2 if cr 7 ^ id, from Theorem 18.41 Lemma [8.51 g(P) = 0 and the estimate 
|V((t)| < Eg =0 (a,u = 0) we immediately obtain: 
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Proposition 8.6 Assuming (17.251) with n < 12 ^ 9d we have 


Y 1^)1 < O^aH 4 * 3 ^" -0 ^) (8.13) 

CTes fc 

a^id, 

for A <C 1. □ 

Proof of Lemma \8.5[ Notice that 1(a) = k only if a = id, for all other permutations 
£(a) < k — 2. We shall prove that, for any £, 

#{a 6 & k : £((t) =£}< 2(2k) k ~ l . (8.14) 

Then (|8.12p follows by recalling k — £(a) = deg(cr) and k < I\ = 0(X~ K ~ S ) and by 
summing up the geometric series 

k 

A 7 deg(,T) < 2 [2A (7_K 

ere©/, m—D 

deg(cr)>D 

To prove (18.141) . let a be a permutation with m ladders of size b\,... ,b m such that 
yT bj = l and bj > 1. If we remove these ladder indices, we have k — £ indices in 
{1,2,..., A:} \ Tg = {ii, * 2 , • • •, ik-(}- The permutation a induces a unique permutation 
a* E &k-i on the indices of this set by a*(j) < a*(j') iff a(ij) < a(ij/). 

Let I* := {0, i\,i 2 , ■ ■ ■ ,ik-e} and set u := |/*| = k — £ + 1. Clearly, the top of any 
ladder belongs to I* and each element of I* can be the top of at most one ladder. We 
assign the length bj of the ladder to its top and for simplicity, we assign the value zero 
to any other element of I*. Thus we obtain numbers b\,..., b u with bj = £ an d 

bj >0. 

If the permutation cr* and the numbers b\,...,b v are given, then we have 2 V ways 
to reconstruct the original permutation a. To see this, first notice that within I* the 
ladder-tops are identified by the condition bj > 0 and the corresponding ladder in a 
can emanate either “to the right” or “to the left” down from its top on the graph of 
a. Once this choice is made, the permutation a can be uniquely reconstructed from a* 
by inserting ladders of given length starting from their tops. Therefore, the number of 
permutations cr E &k with £ ladder indices is bounded by 

V 

2 v (k-£)\x#[(b 1 ,b 2 ,...,K) : Y b j=t,bj> o} 

j = 1 

< 2 V (k - £)\^ < 2(2 k) k ~ e . 

This completes the proof of Lemma 18.51 □ 
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9 Estimates on Feynman graphs with nontrivial 
lumps 

In this section we estimate V (A, a) for a general partition A. We start with a definition. 

Definition 9.1 (%) Let A E Ak- Set a v := \A U \, v E 7(A), to be the size of the v-th 
lump. Let 

S( A):= U ^4- 

*ei( A) 

ci j/ ^ 2 

be the union of nontrivial lumps. The cardinality of this set, s( A) := |/S'(A)| , is called 

the degree of the partition A. 

(ii) Let A E Afc and a E (5*.. The number 

q{ A, er) := max {deg(cr), -a(A)| (9.1) 

is called the joint degree of the pair ( a, A) of the permutation a and partition A. 

The goal is the following generalization of Proposition 18.61 that includes summations 
over non-trivial lumps. 

Proposition 9.2 We assume (17.251) . Let D > 0 and s > 2 be given integers, let 
q := max{D, ^s}. For any k < 33t ^ 36 we have 

Z(k,D,s):= J2 E sup£; (% (A,a, U )|c(A)|<C'(A§-( 6 +^) K - 0 W)VogA| 2 . 

<ree k AeA k u >9—^ 

deg(<r)>£l s(A)>s 

(9.2) 

We recall that for the continuum model non-trivial lumps do not contribute since 
c(A) = 0 unless each A u is trivial, \A V \ = 1. Thus (19.21) can be proved directly from 
Theorem 18.41 Lemma 18.51 and Lemma 19.51 below with a somewhat better exponent. We 
will give a proof for the general case that is valid for both the discrete and continuum 
models. 

Proof of Proposition \9.A The following lemma shows that any even partition P E Vk 
can be generated by a permutation with high degree, depending on the size of nontrivial 
lumps. The proof will be given later at the end of this section. 

Lemma 9.3 For any even partition P E Vk there exists a compatible permutation a = 
<t(P) E &k{ P) such that 

deg (?) > ^s(vr(P)) . (9.3) 
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Corollary 9.4 Given a £ &k and A £ Ak, we have, for k < e+33d 

sup E^ g (A,a, u) < C|logA| 2 (A5-( 1+ ^ d ) K -°(' 5 )) 9(A ’ <7) . (9.4) 

Proof of Corollary \9.f\ We define a permutation a* := a* (A, a) as a* := a if 
deg(cr) > ^s(A), and a* := <t(P(A,ct)) otherwise. By Lemma [iOl we have deg(cr*) = 
g(A,cr). Clearly P(A, a) = P(A,cr*), in particular, E^ g (A, a, u) = E^ g (A, a*, u). 

We wish to estimate the value of an arbitrary partition A by that of the trivial 
partition Ao- We can artificially break up lumps and using the auxiliary momenta to 
account for the additional Kirchoff rules. We describe this procedure in full generality 
for any circle graph. We will call it Operation I because further similar operations will 
be introduced in the companion paper m- 

Operation I: Breaking up lumps 

Consider a circle graph on N vertices (Section 17.11) . Given a partition of the set 
V \ {0,0*}, P = {Pp : p £ J(P)} £ Py, we define a new partition P* by breaking up 
one of the lumps into two smaller nonempty lumps. Let P u = P u iOP v ii with P u tP\P u » = 0 
and P* = {P u t , Pu", P fl : p £ /(P) \ {y}}. In particular /(P*) = /(P) U {v', v"} \ {p} 
and m(P*) = m(P) + 1. 

Lemma 9.5 With the notation above, we have 

P Wg (P, u ,a)< [ drfi W9 (P*,u*(r,y),a) , 

J|r|<ATC 

where the new set of momenta u* = u *(r,v) is given by u* := u fl , p £ 7(P) \ {p} and 
u*, = u v — r, u* v ,, = r. In our estimates we will always have N < 2 K and then 

sup-E(* )ff (P,u,a) < A sup (P*, u, a) 

U U 

with A := [CK(] d = 0(X~ 2dK -° W) (see (S2D and 

Proof of Lemma \9.5l The break-up of the lump P u corresponds to 

d( ^ ±w e -u2)= dr 5^ ^ ±w e - u u + r^)8 ( ^ ±w e -r s j. (9.5) 

eeL±(P,) J e£L±(P v i) eG X±(P„„) 

Note that L(P U ) C L(P U >) U L(P U «) and for any edge e £ ( L(P U /) U L{P u n)) \ L(P U ), we 
inserted an extra w e — w e in the left hand side of (19.51) . Note that the property m on 
the sum of auxiliary momenta is preserved. The integration in (19.51) can be restricted to 
M < IV£ since \w e \ < ( for all e. □ 
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We return to the proof of Corollary 19.41 We apply the break-up operation until all 
lumps become trivial (singlets). This requires not more than s(A) steps, and by using 
Lemma 19.51 we obtain 

sup £ , (*) 5 (A, a*, u) < A s(A) sup£(*) 5 (A 0 ,cr*,u) 

U U 

with A := [CK(f] d . Then (19.411 immediately follows from s(A) < 2q(A,cr) and from 
Theorem 18.41 This completes the proof of Corollary 19.41 □ 

Proof of Lemma \9.HA The inequality (19.31) is trivial if the partition has no nontrivial 
lump, i.e. s( 7 r(P)) = 0, so we can assume s( 7 r(p)) / 0. For any a G 6 *. we define the 

set of internal ladder indices as 

It = IK?) '= ■ |A _1 (* - 1 ) - d _ 1 (z)| = | a~ l {i + 1 ) - A _ 1 (i)| = 1 } 

where a is the extension of a as before. The indices i — 1, i +1 are called the protectors 
of the internal ladder index i G If . They ensure that the index i is neither the bottom 
nor the top index of the ladder. 

We first claim that for any er £ &k, cr 7 ^ id, we have 

k-\If(a)\ <2deg(a) (9.6) 

To see this inequality, we first recall the definition of a ladder and its bottom and top 
indices from Definition 18.31 Since every ladder has a unique bottom and top index, 
that are not internal ladder indices, we see that the sets If, lb and It are disjoint 
subsets of {0,1,..., k + 1}. Since I v ,If C {1,2,..., k}, I v n If = I v D It = 0 and 
k + 1 £ It, we obtain | If\ + \I t \ + \Ib\ + \I V \ h\ < k + 1 + 7, where 7 is the characteristic 
function of the event k + 1 £ h , i.e. 7 := \{k + 1} n h \- Moreover, \I t \ = \ h \, thus 
\If\ + 2\Ib\ + \I V \ h\ < k + 1 + 7. Notice that 1% C If U (If, \ I v ), and k + 1 0 h,I v , 
therefore t < \If\ + \Ib \ I v \ — 7 = \If \ + \h\ — \Iv C h\ — 7- From the last two inequalities 
and from deg(cr) = k — I, we obtain 

k - \If\ < 2 deg(cr) + 1 — 7 — \I V \ I b \ — 2\I V C I b \ < 2 deg(<r) + 1 - \I V \ . 

Since a 7 ^ id, we have |/„| > 1 and which proves (19.61) . 

Next we will show that there exists a compatible permutation a G ©^(P), 

i?(a)na(5(A))=0, ( 9 . 7 ) 

where we set A := 7 r(P) for simplicity. Since 5(A) / 0, If (a) 7 ^ {1, 2,... , 7’}, i.e. 
<7 7 ^ id. Combining (19.61) with (19.71) and with the fact that both If and <r(5(A)) are 
subsets of {1, 2,... , k}, we obtain (19.31) . 
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To construct a satisfying (19.71) . we apply a greedy algorithm. Since P is even, ©^(P) 
is nonempty and we pick a do £ Sfc(P). If (19.71) is not satisfied for <to, then some internal 
ladder index i is in the image of a nontrivial lump A £ A; i = cro(i'), i! £ A. Let j’ £ A 
be another element of this lump. Flip the permutation do on these two elements, i.e. 
define a new permutation a i such that oq [i!) := do (j') = j , <7i(/) '■= do(i / ) = i, and 
(T\(r) := do(r) for any r / i,j. Clearly o\ £ ©fc(P), di(s(A)) = do(s(A)). We claim 
that 

|/;(di) n di(5(A))| < |/;(d 0 ) n do(5(A))| , (9.8) 

i.e. the total number of internal ladder indices in nontrivial lumps decreased. Continuing 
this flipping process for di etc., we obtain a permutation a satisfying (19.71) . 

To see (|9.8D we note that after the flip i is not an internal ladder index any more. This 
is clear if j ^ i — 1, i + 1; in that case the points (d _1 (i — 1), i — 1) and (d _1 (* + 1), i + 1) 
have not changed and they would uniquely fix the location of an internal ladder index in 
between. The preimage of the index i has moved out from this position, dj” 1 (i) A d ( ^ 1 (z). 
The index j however would not become internal ladder since d^ 1 ^') = i' is between 
dj _1 (i — 1) and d^(i + 1), but j is not between i — 1 and i + 1. It is easy to see that 
the fixed points (d^ 1 (i — 1), i — 1) and (d _1 (i + 1), i + 1) also prevent any other indices 
from becoming an internal ladder index after the flip. This could only be possible if due 
to the new point (cfj" 1 ^'), j) = (d ( ^ 1 (z),j), one of the neighbors of j, say j + 1, would 
become an internal ladder index. It is easy to see that then j + 1 must be equal to i — 1 
and the other protector of the new internal ladder index j + 1 must be i. In this case 
i — 1 was already an internal ladder index before the flip as well, so no new internal 
ladder was created. 

A similar but simpler argument shows that if j = i — 1 or j = i + 1, the number of 
internal ladder indices also decreases. This completes the proof of Lemma 19.31 □ 

We continue the proof of Proposition 19.21 Given a £ &k and A £ Ak, we recall 
the permutation a* = cr*(A, a) defined in the proof of Corollary 19.41 We also note that 
s(A) < 2deg(cr*). Hence 

E{k,D,s)< l c ( A )l J2 J2 sup E^ g (A, a, u) . 

s(A)>s deg(er* )>q cr* (A,cr)=<r* 

Note that a £ &k{ P(A, cr*)), so by (|7.20l) the summation over a contributes by a factor 
of at most a u \ and we obtain 

Z(k,D,s)< Y ( II a “"~V ! ) sup£ ( * )!? (A,(7*,u) . 

<r*ee k AeA k vel(A) U 

deg(a*)>q 2<s(A)<2deg(a*) 


48 




We also used the estimate (16.111) . By using (19.41) and deg(<7*) = q(A,a), we obtain 


E(k,D,s)< •£ (CAJ-< 1 +T«-om) d0s( "‘ ) y ( JJ a ;~- 2 aj). 

A eA k veI{A) 

deg(cr*)>q s<s( A) <2deg(cr* ) 

We introduce the notation 

* * 

^fM:= y fM ’ n/M- n fM- 

vei(A) vei(A) 

a u >2 

First we £x the sizes of the nontrivial lumps a v > 2. Given these sizes, the number 
of A partitions is bounded by 

f k\fk — aA fk — ai — a,2\ ^ k\ k^- 

V«1/ V «2 / V a 3 ) " ' ~ {k - Y,* av) ] -W* a v \ ~ W* a v \ ' 

Recalling s(A) = Yl* a vi and that s(A) < 2deg(cr*), we have 

* 

S(k,D,s)< Y (Ck 2 \ T ) Aeg( ' a *' > Y (n«r 2 ) (9.9) 

^*ee k a „ : J2* a u < 2deg((T*) 

deg(u*)>ij 

with r < i — (1 + ^d)n. We use the bound a a ~ 2 < C a ~ 1 (a — 1)!. To estimate the 
summation over a u ’s we use the following inequality. For any fixed m, H we have 

# * 

£(HK-!)!)< (*-!)! (9.10) 

where the summation # is over all sequences (ai, 02 , •.., a m ) of positive integers at least 
2, whose sum is H. The proof of (19.101) is easily obtained by induction on rn from 


H -2 h -2 

J>-l)!(ff-a-l)!< £(ff-2)!<(tf-l)!. 

a=2 a=2 

Summing (19.101) over all H < 2deg(cr*) and m < H/ 2, we obtain the bound 

y (n a “"' 2 ) ^ 2 t 2de g(o]! < ( ck ) 2deg (-*) 

Ok : a„<2deg(o-*) 

for the 0 ^ summation in (19.91) since deg(i7*) < k + 1 by definition. 
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In summary, we obtain from (19.91) 

E(k,D,s)< ^ (°k 4 X T ) dee{a) 

CTes fe 

deg(<r)>ij 

Recalling that k = 0(A _K_<5 ), we can apply Lemma 18.51 with 7 = r — 4 (n + S) as long 
as 7 > k + 5. For sufficiently small positive <5 this gives the condition k < in 

Proposition 19.21 and the estimate (19.21) . □ 


10 


Proof of Theorem 


[833 


For any (k + 1) x (k + 1) matrix M we set 

r Y r 


E(M ) :=supA 

U 

k+1 


2k 


dad/3 / d//(p) |R(p)| \B(Mp + u)||V’(Pi)IIV’(Pi +«i)| 


-Y 


X 


n 


1 


1 


\a - u{pj) - irj | |(3 - ui([M p + u]^) + irj\ ' 
The key step in the proof of Theorem 18.41 is the following Lemma: 


( 10 . 1 ) 


Lemma 10.1 Suppose (|7.25l) . assume that k < and 6 < 5 (k) is sufficiently small. 
Then for any a £ 6^ we have 


E{M[a)) < C'^A^ _ ( 1+ ^ <i ) K_0( '' 5) ^ deS(CT) | log 77 I 2 , ( 10 . 2 ) 

where the matrix M(a) was defined in (18.71) and the degree of the permutation, deg(cr), 
was defined in »• 


From (18.41) . (18.81) . clearly 

e 2tri 

sup E g= o (a, u) < ——2 E(M(a)) (10.3) 

u (27T j 

after integrating out all pj variables in (18.41) and by using pi = pi — The estimate 
(jl0.2l) will then complete the proof of Theorem 18.41 for g = 0. 

The proof of Theorem 18.41 for other (but finitely many) g values follows exactly in the 
same way. This requires to slightly redefine B( p) (see (16.51) 1 in the definition of E(M ) 
by allowing the factor {pj+\ —Pj)~ 2d instead of B(pj + \ —pj) at a few places exactly as in 
the definition of Mg (see ()7.3D ). As we will see along the proof of ( 1 10 . 2 1 ) , this change will 
require using the less precise bound (|3.13j) with a = 0 and h(p—q) = (p—q)~ 2d instead of 
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the more accurate estimate (13.151) at most g times. Each time we lose a constant factor 
compared with the proof for g = 0. Since g < 8, this results only in a constant factor. 
Finally, the proof for the truncated E -values requires to define a truncated version of 
E(M), where the last product in (110.11) runs only up to j = k, i.e. the last a and /3 
denominators are not present. It will be clear from the proof of Lemma 110.11 that the 
same bound holds for the truncated version of E(M) as well. □ 


10.1 Pedagogical detour 

The size of the multiple integral in (110.11) heavily depends on the structure of M = 
M(a). Before we go into the algorithm to evaluate this multiple integral, we present two 
calculations, that introduce the techniques that we are going to use in the actual proof. 
The second calculation also provides the bound (110.21) . hence (18.111) . for the case of the 
trivial permutation deg(cr) = 0. 


10.1.1 Method I. Pointwise bound 


The most straightforward bound on (110.11) estimates all but one of the /3-denominators 
by L°° norm 


Ap %j I P ~ u([Mp + u]j) + ip | 


< 7 / 


-1 


(10.4) 


It would be possible to estimate this denominator by C A~ 2 apart from a neighborhood 
of zero using (13.71) and treat the |p| ~ 0 regime separately. Here we choose the simplest 
argument and we do not optimize for the best possible exponent k. 

We use (110.41) k times to obtain 


E(M) < 


C X 2 \k 


1 


sup 

U , 


d]i{j) k+x 


y 


(Pk+ 1 + u k+ i) 2d J J—Y |« - u{p k+ 1 ) 


/.../d PI ...d M n 


\B(Pj+i ~Pj)\ 
L 1 ! |a - w(pj) - iij\ 


dad/3 

ir/ 11/3 - u}(p k+ i + Uk+i) + ir]\ 


Here we used the following bound for any q = (qi,, q k+ 1 ): 

|S(q)ll^o(9i)| <C k (q k+1 )~ 2d (10.5) 


to obtain the decay in p k +i = p k +i + u k +i- We integrate out pi,P 2 , ■ ■ ■ ,p k using (13.121) . 
then we perform the da, d/3 integrals and finally dp k+ \ to obtain 

E(M) < (C\ logA|) fc+2 (A 2 r j ~ 1 ) k . (10.6) 

The estimate (110.61) is off by a factor (A 2 r/ _1 ) fc = ( X~ K ) k because we did not use the 
stronger estimate mentioned after (110.41) . We also collected many logarithmic factors 
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and the constant is not optimal. We note that in the typical term k ~ X 2 t ~ X~ K 1, 
so even an error C k may not be affordable. To improve this estimate, for a typical matrix 
M, we will not use the pointwise bound (110.41) for all /3-denominators. We will carefully 
select those /3-denominators whose singularities cannot overlap with other singularities, 
hence they can be integrated out at a | log 77 1 expense instead of r/ -1 . 

Before we explain this algorithm, we show another method to estimate E(M). It 
practically estimates E(M) by E(I), i.e. by the ladder graph, that can be computed 
more precisely. The same calculation will be important when evaluating embedded 
ladder graphs. 


10.1.2 Method II. Successive integration scheme for ladder graphs 

We separate all but one a and /3 denominator by a Schwarz inequality. We obtain 

dad/3 


E(M) < X 2k sup / dp(p) 


Y 


—Y l« - w(pi) - irj 11/3 - u(pi + hi) + ifj\ 


fc+i 


\Mpi)\ 2 n 


3=2 


\B{pj ~Pj~i)\ 2 
| a — Zo(pj) — ip\ 2 


k -\-1 


+i^o(9i)i 2 n 


3= 2 


\Bjqj ~ qj- 1 )\ 2 
\/3 - u(qj) + ir/\ 2 


with the shorthand notation q := Mp + u. Because M is an invertible matrix with 
determinant ±1 (Proposition 18.21) . the contributions of the two terms in the square 
bracket are identical up to exchange of a and (3. To estimate the first term, we use 
iteratively (13.151) and (13.131) (with a = 1/2) to integrate out pk+i,Pk, ■ ■ ■ 1 P 2 (in this 
order): 


A 2 


\B(pk+i ~Pk )| 2 
|a - TD{p k+ i) - ip | 2 


dpk+i < [l + CoA 12k (A + |« —w(p fc )| 1/2 ) 


(10.7) 


2 r 


p^=7 P s k 1 ^| 2 [ 1 + (:7 oA 12k (A + |a - u>{p k )\ 1/2 )\dp k 
|a - u(jpk) - ip\ 2 L J 

< (1 + CA 1_12k ) [l + C 0 A- 12k (A + |a - c 4Pk-i)\ 1/2 ) 


( 10 . 8 ) 


etc, with C := C 0 ( 1 + C' 1/2 11-S 2 1| 2 rf,o)- In the last step we use only (13.131) once for a = 0 
and once for a = 1/2: 


A : 


/ 


\B(p 2 -Pi)| 5 


l + C 0 A- 12K (A + |a-ta(p 2 )| 1/2 ) 


|a - u(p 2 ) - ip | 2 
Then we integrate da d/3 and finally dpi to obtain 

E(M) < C{ 1 + c , a 1 ” 12k ) A: | log A| 2 < C\ log A| 2 


dp2 < C . 


(10.9) 


( 10 . 10 ) 
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by using k < K <C A -1 " 1 " 12 ” as k < 1/13. 

We note that this method also gives a robust bound for the truncated -E-value, since 
the truncation means that Lemma 13.31 is used only A; — 1 times. Summarizing, we have 
proved 

Lemma 10.2 We assume (17.251) and k < 1/13. Then 

sup sup.E(cr, u) < d7| log A| 2 

o-eSj, u 

sup sup E*(a, u) < C7A 2 1 log A| 2 . □ 

0-e6 fe u 

10.2 Choice of the integration variables 

Before we start the proof of Lemma riO.il we explain the main idea. We use a combination 
of Methods I and II. We will assume in the sequel that a ^ id. The lemma for the trivial 
a = id case has been proven in (110.101) . 

Note that each factor in the integrand in (110.11) is almost singular on a set of codimen¬ 
sion one of the form {a = Re u(pj)} or {(5 = Re ui([M p + u]j)} in the high dimensional 
space of integration, (R rf ) fc+1 x R 2 . The singularities are regularized by ?y and Im u 
and the two regularizations always have the same sign. The matrix M may enhance 
the strength of these singularities by forcing these “almost singularity” sets to overlap. 
For example, in the ladder diagram with u = 0, we have p = p, hence the singularity 
sets are pairwise identical if a = (3. Therefore singularities of quadratic type necessarily 
occur. It is expected that this is the only mechanism that creates relevant overlaps 
of singularities. Hence, ideally, one would integrate out the ladder momenta with the 
precise bound (13.151) . This would remove all denominators with indices j £ Ig and the 
remaining integral should be of 0( A 2 ( fc ~^) with possible logarithmic corrections, i.e. one 
should gain a factor A 2 from each non-ladder index. 

However, the singularity sets of the remaining denominators may still intersect on 
higher codimensional manifolds. For example, since 

r _ i _ i c 

S aS j la - u(p + q) - iri\ \(3 - w(p) + irj\ ~ \q\ + rj' 

the integration of these two propagators develops an almost-singularity at q = 0. We will 
call such factors point-singularities, although they are regularized at a very short scale 
ry. After integration, one point singularity is harmless as (|g| + ry) _1 £ L] oc (dq), how¬ 
ever, several point singularities may accumulate along the procedure whose simultaneous 
integration may lead to further divergences. 

These enhancements of singularities are expected to have contributions of lower or¬ 
der, but their estimate is not easy. Note that every integration variable pj may appear 


( 10 . 11 ) 

( 10 . 12 ) 
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in many denominators. This interdependence renders the effective /^-estimate of each 
integral practically impossible. A Schwarz inequality (Method II) can remove all cor¬ 
relations between denominators, but the resulting T 2 -estimate is of the same order as 
the main (ladder) term and we would not gain anything from the higher degree of the 
permutation a. 

The idea is to estimate many, but not all /3-denominators in (110.11) in the trivial 
way (110.41) . These denominators are chosen in such a way that the remaining ones can 
be successively integrated out without ever computing an integrand with more than 
two propagators and without collecting more than at most one point-singularity factor. 
This choice of integration variables are determined from M = M(a) by an algorithmic 
procedure that we describe now. 

We first recall that M is a tower matrix, and we define b(j) and t(j) to be the bottom 
and the top of the tower in the j-th column, i.e. 

Hj) := max{z : M tJ / 0}, t(j) := min{/z : M tJ / 0} . 

We also recall the concepts from Definition 18.31 Any non-peak index, i £ I \ I p , clearly 
has the property that i is the bottom of the tower of some column of M, i.e. there exists 
j such that b(j ) = i. This column index can be either j = d' _1 (z) or j = cf -1 (i) + 1 
or both; ambiguity occurs only if i £ is a valley index. For any ladder or slope or 
last index, i £ Ig U I s U Ii as t-> we therefore define c(i ) to be the unique column index j 
such that b(j) = i. For i £ Ii as t, be. i = k + 1, we always have c(k + 1) = k + 1. If 
i £ I v , we define c{i) to be index of the column whose tower is higher, i.e. c(i) := cr~ l {i) 
if t(a~ l (i)) < t(a~ 1 (i) + 1) and c(i) := <7 _1 (z) + 1 otherwise. The other column index 
will denoted by c(i); this concept is defined only for i £ I v . The i £ I \ I p , 

elements of the matrix M are called pivot elements; these will be used for determining 
the integration variables. The c(i)-column will be called the pivot column of i and 
the corresponding tower in this column is called the pivot tower of i. The definition 
guarantees that the pivot tower of an index i £ I s U I v has at least length two. 

Definition 10.3 Let h\ < 1%2 < ■ ■ ■ denote the elements of I v U I s in increasing order. 
A slope index h p £ I s , is called covered if t(c(h p+ 1 )) < hi.e. if the top of the tower 
of the next pivot element of a valley or a slope is higher than h p . The remaining slope 
indices are called uncovered. The set of covered and uncovered slope indices will be 
denoted by I cs and I us , respectively. 

Fig. [6] shows three examples. In the first two cases, the index h p is covered by slope 
and by a valley. The last case is an uncovered slope index. The shaded boxes indicate 
the pivot elements. 

Apart from the ladder momenta, p c (j), i £ Ig, we will integrate out the variables 
Pc(hi)iPc(h 2 )) ■ ■ •> i n this order. We will show that at every integration step only one 
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Figure 6: Two examples for a covered slope index and an uncovered index 

a-denominator, namely \a — Co(p c ( h ^) — ip | _1 , and one /3-denominator, namely \/3 — 
o;(p/j M ) + will contain the integration variable p c (h ), so we will never have to 

estimate integrals with more than two propagators. 

The significance of the covered slope indices is that their integration step yields a 
point singularity that will be integrated out immediately at the next non-ladder inte¬ 
gration step. Therefore covered slope indices do not give rise to accumulation of point 
singularities. 

The covered slope indices, together with the valley indices and the last index will be 
used to gain a \ const factor in the integration procedure. They will not accumulate point 
singularities. The uncovered slope indices and the peak indices will be treated by the 
trivial L°° bound (110.4|) and no gain is obtained from them. That is, we prefer giving up 
the potential gain of a \ const factor from the uncovered slope indices to dealing with the 
accumulated point singularities. Finally, the ladder indices will be treated by Method If 
with the Schwarz inequality. 

Lemma 10.4 The number of uncovered slope indices is at most v, where v = \I V \ is the 
number of valleys. 

This lemma ensures that treating the uncovered slope indices by the trivial bound is 
affordable (at the expense of a worse k), since each valley will result in a gain of a \ const 
factor. 
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Proof of Lemma 10-4\ First we state two claims from which we easily deduce the 
Lemma, then we prove these claims. Recall that the sequence h\ < hi < ... lists all 
elements of I v U 


Claim 1. Let h v E I s be the first slope index in the sequence hi < h ,2 < ■ ■ ■■ Then 
there is a peak index h such that h < h v . 

Claim 2. Let L /t E I us be an uncovered slope index and suppose that it is not the 
last entry in the sequence h\ < h -2 < ..., i.e. there is a next element, h p+ 1 E I s U I v . 
Then there is an intermediate index, h p < h < h^+i, that is a peak index. 

By Claim 2, one can assign a distinct peak index to every uncovered slope index 
h, t (but the last one), namely the next peak index. Moreover, by Claim 1, there is an 
additional peak index smaller than any uncovered slope index. Therefore the number of 
uncovered slope indices is at most the number of peaks. The statement the lemma then 
follows, since the number of peaks and valleys is the same, p = v. 

Proof of Claim 1. Suppose that the statement is false. First note that h u >2 since 
index 1 cannot be slope. Then all indices 1,2,... h u — 1 must be either ladders or valleys. 
But for any valley index j there is obviously a peak index j' < j, since the graph of the 
permutation goes down from (0,0) to (k + 1, k + 1), so following the two branches of the 
graph emanating upwards starting from a valley at least one of them must run into a 
peak that lies higher. Therefore all indices 1,2,... h v — 1 must be ladders and they form 
a ladder, with the index 0 being the top of this ladder. In particular t(c(h )) = h — 1, 
h = 1,2,... h u — 1 (recall that t(c(h)) denotes the index of the top row in the pivot tower 
of h ), so every index 0,1, 2,... h u — 2 is the top of a column of length one. Moreover, 
none of these indices can be the top of any other column, since an index which is the 
top of two different columns must be a peak index. But in that case there is no room for 
t(c(h u )), the top index of the pivot column of h v : since h u was a slope, t(c(h u )) < h v — 2. 

Proof of Claim 2. To see the claim, suppose that none of the indices /i, strictly 
between h p and h p+ 1 , is a peak. These indices cannot be slopes or valleys either since 
hfj, and h p+ \ are two consecutive elements of I s U I v and h cannot be the last index, 
h 0 Iiast j since h < h p+ \ < k + 1. Therefore all indices h with h p < h < h p+ 1 are 
ladders. 

Now we look at h* = t{c{h^ + 1 )), i.e. the index of the top row in the pivot tower of 
h ^ + 1 (see the last picture on Fig. [6]). Since h M was uncovered, < h*. Note that the 
pivot tower of h^ + \ E I s U I v has length at least two, therefore h* < h^ + 1 , and thus 
h* E Ii. Since h* was also the top of a tower with length at least two, h* — 1 must be a 
peak index since it is a row index right above the top of two towers (namely the towers 
in the columns c(h*) and c(/r /i+ i)). Clearly < h* — 1 < h^ + 1 , but / h* — 1 since 
h fi E I s and h* — 1 E I p . Thus h* — 1 is a peak index strictly between h p and h^+i, 
which is a contradiction. □ 
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10.3 Integration procedure 

Our goal is to estimate (110.11) for M = M(a) when a ^ id. We start with defining 


g||| := ?? + min{|g|, 1 }, q £ R d . 


(10.13) 


This is not a norm, but it satisfies the triangle inequality, |||p + g||| < |||p||| + |||g|||. For 
any index set I' C / = {1, 2,... , k + 1} we define the function Uj> as the product of 
those potential terms, \B\, in (110.11) that depend only on momenta {pj : j £ I'}. More 
precisely, 


Uv{ P,u) 


B(Pj +1 ~ Pj) + u\ j+ i - [Mp + u]j) , 

3 3 


where the star indicates a product on a restricted index set. The first product is taken 
over all indices j for which j,j + 1 £ I'. The second product is taken over those j's, for 
which Mj+ib = Mjb for all b 0 V. 

For any \I'\ x (k + 1) matrix M and any vector b = (&i, 62 , • ■ ■, &fc+i) £ R fc+1 we define 

E(l',M,b ) :=A 2fc sup [ [ dad/3 sup [ IT dp(pj) Ujr(p,u)-rr ——--- 

U,D J J-Y Pj : jgr J j eI , It ■ p + V\l 

x f rr _ 1 _I_ 

\j e j, l« - Pj ) - iy\ \P - w([MP + u]j) + iri\ 

Here v £ is an additional dummy momentum and b ■ p is defined analogously to (18.51) 
as 

b • p := b\p x + b 2 p2 + • • • + bk+iPk+i ■ 

We will also use the notation E(I', M, 0) defined exactly by the same formula as (I1U.141) 
without the factor |||b • p + u ||| _1 in the integrand and without the supremum over v. 
With a slight abuse of notation we will refer to this case as chosing the “empty vector” 
b = 0. 

For each h £ I we will define an index set I^ h \ a matrix M^ h \ and a vector b^ 
such that E(I^ h \ M^ h \ b^) is the remaining integral after the h-th integration step. 
Each step consists of eliminating one /3-denominator; in the h -th step we will eliminate 
|/3 — w([Mp + u]/j) + iql^ 1 and sometimes one or two additional a-denominators. The a- 
denominators will always be eliminated by integrating out the corresponding p?-variable. 
The order of integrations will be such that the integrand in E(I^ h \ M^ h \ b^)) does not 
depend on any pj 0 I^ h \ In particular, the supremum over all pj 0 1^ in the definition 
HO.141) will be superfluous in the integrals appearing in our iteration. 


(10.14) 
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The elimination will be done differently, depending on the type of the index h. If 
h € I p , then the trivial L°° bound (110.41) will be used. If h is an uncovered slope index, 
then again the trivial bound will be used, and, additionally, the variable p c (h) will be 
integrated out. If h is a slope or ladder index, then variable p c (h) will be integrated out. 
If h is a valley then we integrate out both variables, p c ^ and Pz(h)- Finally, in the last 
step, h = k+ 1, we will integrate out the j3 variable. We will bookkeep this integration by 
removing the corresponding h rows from M and we will say in short that these rows have 
been integrated out. At the end we will have to integrate out the remaining p^-variables, 
but then only one ce-denominator will be present so the integration can be done easily. 

Point singularities are created by integrating out a covered slope index, but there 
is never more than one singularity present. An existing point singularity will not be 
touched by the ladder integrations and it will be changed to another point singularity 
when a covered slope index is integrated out. In all other integration steps an existing 
point singularity will be removed. 

Now we define M^ h \ and b^. Let := M = M(a), 1^ := I = {1,2,, k + 
1} and = 0. In particular, note that 

E(M) < H^oIlL £(/ ( 0 ) ,M(°\bW) . (10.15) 

For any h > 1, define to be the matrix M = M(a) with the first h rows, with 

indices i = 1, 2,... , h, removed. For h > 1, let 


I {h) :=I\ ^{c(h') : h' <h,ti € I \I P } U {c{ti) : ti < h, ti G 


be the index set of columns that have not been integrated out up to the h -th step. 

The vectors b 1 ^ = {J}[ \ ..., are defined as follows. For any given h G I, let hp 
be the largest element in the listing h\ < h ,2 < ■ ■ ■ of I v U I s (see Definition 110.31) such 
that hp < h. If there is no such hp or hp 0 I cs , then we set b^ = 0. If hp € / cs , then 


(h) f M h/x;j if j ± c(h) 

J ' 1 0 if j = c{h). 


(10.16) 


We collect a few information about b^d that directly follow from the definition and 
from the structure of M. 

(i) Suppose that ^ 0, and hp as above (i.e. hp < h and it is the largest element 
with this property). Since hp e I cs , i.e., t(c(hp+ 1 )) < hp, we therefore have 

b%l +i) = ±i, hp<h<hp + 1 , hp G I cs . (10.17) 
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(ii) For any fixed p, b^ is constant as h runs through the interval {h p ,..., h p+ \ — 1}. 

(iii) Consider the pivot column c(i) of a ladder index i e Iy. Since the only non-zero 
entry in the c(i) column of M is M l c ^. we have that My ^ = 0 for all p. But 
then we see that the entries corresponding to any ladder pivots are zero in the 
vectors: 

bc(i) = 0 f° r any i £ Iy and any h such that b^) / 0. (10.18) 

To describe the integration, we explain how each case, depending on the index type 
of h, is estimated. We also introduce the short notation 

E(h) := E(I {h \M w ,b {h) ) , h = 0, 1,... ,k + l. 


Case 1: h G I p . In this case no integration is done, = I^ h 1 ( and b^ = b^ h . 
We estimate the denominator | (3 — u([Mp + u]h) + ip\~ l trivially as in (|10.4j) and we 
obtain 

E(h — 1) < rf~ 1 E(h) , h € I p . (10.19) 


Case 2: h € Ig. Let h, h + 1,... ,h + t — 1 £ Iy be & ladder, i.e. a maximal sequence 
of consecutive ladder indices, i.e. h + r 0 Iy with some r > 1. We will integrate 
out the whole ladder in this step, so we can assume that h — 1 is the top of a ladder, 
h — 1 6 It- By definition of the ladder, the pivot indices, c(/i), c(h + 1),..., c(h + r — 1) 
are consecutive numbers, which, for definiteness, are assumed to be increasing (the other 
case is similar). Let c = c(h). Because h ,..., h + T — 1 are all consecutive ladder indices, 
[Mp]j = Pi + w for all c < i < c + r — 1 , with a vector w = w(j>) that does not depend 
on {pj : c<j<c+r— 1}. Note that b remains unchanged as k runs through the 
ladder indices: b ^ = b for all k = h, h + 1,..., h + r — 1. We claim 

E(h 1) < C\~ 2r E{h + t - 1) , h-l e I t (10.20) 


with a constant C independent of r. This inequality entails integrating out the variables 
p c ,p c +i ,... ,p c+T -i (in this order) to remove the rows h,h + l,... h + T — 1. This requires 
estimating the following integral 


T := sup 

a,(3,w * 


. c+t—1 

n 

j=c 

C+T — l 


dp(pj 


a - u){pj) - ip 11/3 — iv(pj + w + uj) + ip\ 

x II l-®(Pi+i _ Pj)\ l^(Pi+i + “i+i - Pi ~ «i)l 


( 10 . 21 ) 


J=C 


because the possible point singularity present in _E(/(^ _1 ), b^ -1 ^) does not con¬ 

tain any of the integration variable by (110.181) . Due to the ladder structure it is easy to 
see that Uph-i) indeed contains all the B factors shown in (110.2ip . 
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After a Schwarz inequality we use an iterative integration scheme similar to (110.71) - 
(I10.9D to estimate the r-fold multiple integral 


f \B(Pc+t Pc+t— i)| ^m(Pc+t— l) 

J | a — aJ(p c+r _i) — ip\ 2 


f \B{p c +2 - Pc+i)\ 2 dfi(p c+1 ) 

J \a -uJ{p c+ i) -ip\ 2 

< [A _2 (l + C'A 1_12k )] t_1 


|-6(p e+r _i - p c+T -2)\ 2 dfi(p c+T - 2 ) 

| a - uJ(p c+r _ 2 ) - ip | 2 

r \B(p c+ i ~ Pc)\ 2 dp(p c ) 

|a - uj(p c ) - ip | 2 


/ 


I B(p c+r — p c _|_ T _ i) | 2 d/r (p c+T _ i) 


|a - w(p c+T _i) - m?| 2 

< C[A" 2 (1 + CA 1 " 12k )] t_1 < CA" 2 ^" 1 ) . 


1 + CqA 12k (A + \a — cu(p c _|_ T _i)| 1,/2 ) 


In the last integration we use (13.131) with a = 0 and a = 1/2 and we also used that 
r < if <C A -(1 ~ 12k) . This gives 

X < C\~ 2t . ( 10 . 22 ) 


Case 3: h € I us . We have i^ = \ (c(/i)} and b^ = 0. We start proceeding 

as in Case 1 by first estimating the denominator |/3 — oj([Mp + u]^) + ip\~ l trivially as 
in (110.41) . Then observe that the variable p c (h) appears only in one propagator, namely 
in |a — cJ(p c (h)) ~ ip | _1 - Moreover, if b^ _1 ) ^ 0, then p c p^ appears also in the point 
singularity, since b^ h] L} = ±1 by (110.171) . 


Before we integrate out p c (h), we estimate all B factors that contain this variable by 
supremum norm. The number of such factors is at most four so they can be bounded 
by a constant factor C. Then, to integrate out p c n^, we either use the bound (with 
P = Pc(h )) 


sup 

|a|<y. 


d pijp) 


| a — uj(p) — ip | 


< CC d 2 \ logr/| 


(10.23) 


if b^" 1 ) = 0 , or we use the bound 


sup sup 
\a\<Y r . 


d pip) 


|a — uj(p) — ip | |||p — r\\ 


< C( d ~ 2 \ log 77 1 


(10.24) 


if b^ _1 ) 7 ^ 0 (recall ( = \~ K ~ 3S from (17.25^ 4 . The bound (110.241) follows immediately 
from (1A.1I) and (IA.5I) . the bound (110.231) is weaker since sup p |||p||| < 1 + p (the bound 
(110.231) can be improved to C\ logry| but we will not need this fact). In summary, we 
have proved that 

E{h - 1) < CC^p- 1 1 log 77 1 E{h) , helus- (10.25) 
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Case f: h G I cs . In this step we will have 1^ = I^ h \ {c(h)} and will be the 

h -th row of M but with the pivot element (h, c(h )) changed to zero, i.e., = 0. Since 

index h is among the listing hi < /12 < • • • of the elements I v U I s (see Definition 110.31) , 
let h = hfj, and let h! := h jl+ ] be the next element of I s U I v . Since h is a covered slope 
index, we have t{c{h')) < h and thus we know that the (h,c(h')) matrix element of M 
is ±1. Therefore the variable p c (h') appears non-trivially in ph = [M p + u]^: 

Ph = \Mp + 1 % = ±p c{ h) ± p c ( h/ ) ±... ( 10 . 26 ) 


After estimating those B factors by supremum norm that contain p c ( h ), we will integrate 
out p c (h) by using the following lemma proved in the Appendix: 

Lemma 10.5 Let A 3 < p < A 2 , and |g| < CA -1 . Then we have 


f d p(p) 1 ^ Cp 1 / 2 ( d 3 

a/,r J \a- ZJ(p) - ir)\ \(3 - u{p + q) + ip | ||| p - r||| “ ||M|| 

Without point singularity we have the following improved bound 

f _ d pip) _ < C( d ~ 3 \logp \ 2 

8 a$ J \a- w(p) ~ ip | \(d - w{p + q)+ ip\ ~ |?| 


log ? 7| 2 


(10.27) 


(10.28) 


The first estimate (110.271) is used if b^ h b fz 0 ; i. e . if there was a nontrivial point 
singularity present. By (110.171) . this nontrivial point singularity must originate from 
the previous covered slope index, in particular /i M _ 1 G I cs and 1 = ±1, i.e. the 
point singularity contains the variable p c (h)- The new point singularity, obtained on 
the right hand side of (110.271) as |||g||| _1 , is exactly the linear combination appearing 
in (110.261) without the integration variable, p c (h)- Therefore the p-dependence of this 
point singularity is given by b^ (up to an overall sign that is irrelevant since ||| ■ ||| is 
symmetric). If b^ -1 ) = 0, then we can use the second estimate (110.281) and otherwise 
we follow the same argument as in the b^ _1 ) 7 ^ 0 case. 

The net effect of the elimination in Case 4 is 

E{h- 1) < Cp- 1/2 ( d ~ 3 \logp\ 2 E(h) , h G I cs . (10.29) 


Case 5: h G I v . In this step we will have \ {c(h),c(h)} and b^ = 0. 

Recall that c(h) and c(h) are the two columns between which the valley is located, 
i.e. {c(h),c(h)} = {d _1 (h ), a~ x (h) + 1}. We integrate out both p c (h) an d Pc(h) after 
estimating the B factors that are involved by their supremum norm. This eliminates 
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two a-denominators and the |/3 — uj([Mp + u]^) + ip\ 1 denominator. It follows from the 
definition of the valley index that 


[M p + u]h = ±p c (h) T Pc(h) ± • • • 


i.e. the /3 denominator contains both integration variables. If b^ h ^ = 0, i.e. no point 
singularity is present, then the prototype of this integral can be estimated as follows: 


sup 

n,m <y 


dp{p)dp(p') 


< C( 2d ~ 5 \ logr/l" 


|a — co(p) — ip | \a — ui(p') — ip\ \(3 — u>(p — p' + q) + ip\ 


assuming |g| < C A 1 . This estimate follows from (110.28jl and (110.241) . 

If a point singularity is present, / 0, we know that 6^^ = ±1 by (110.171) . 

therefore at least one of the integration variables appear in it. Depending on whether 
c(h) also appears or not, we have one of the following two prototype estimates: 


sup sup 

<Y r 


_ dp(p)dp(p') _I_ 

a — l J(p) — ip\ |a — ZJ(p') — ip\ \/3 — ui{p — p' + q) + ip\ |||p — r 


< C7 ? - 1/2 C M_5 |log?y| 3 


(10.30) 


or 


sup sup 

\»m<y r 


dp(p)dn(p') 


\a — u)(p) — ip\ |a — u(p') — ip\ \/3 — uj(p — p' + q) + ip\ |||p ± p' — r|| 


< Cp~ 1 / 2 C 2d ~ 5 \ logp| 3 , (10.31) 

assuming |g| < CA^ 1 . Both inequalities are obtained by first applying (110.271) for the dp 
integral, then applying (|10.24[) for the dp' integral. 

In summary, the effect of the elimination in Case 5 is 

E(h-1) < Cp~ 1/2 C 2d ~ 5 \logp\ 3 E(h) , he I v . (10.32) 


We eliminate the row indices one by one in increasing order according to Cases 1-5. 
The total estimate is 

E( 0) < (r ? - 1 ) i/pl A- 2£ C |/t| (CC d “ 2 r? _1 |logr/|) |7us l(Cr ? - 1 / 2 C d_3 |log? 7 | 2 ) |7cs| 

x(Cp~ 1/2 C 2d ~ 5 \logp\ 3 ) lIvl E{k + 1) , (10.33) 

where we recall that It is the set of indices that are tops of a ladder, in particular, Case 
2 is applied \It\ times. Since the top of the ladder is not a ladder index and not a valley, 
so it is a peak, slope or maybe 0, thus 

|-It | £ |/p| + |/ s | + l = u + s + 1 < 2v + s 
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taking into account that p = v and v > 1 since a ^ id. From (110.331) and from 
|I CS | + 1 1 us | = |/ s | = s it then follows that 

£(0) < C v+s+1 \~ n ri-^ 3v+s+ ^C d{ - 2v+ ^E(k + 1) (10.34) 


(we also used that | logr/| 3 <C ()■ Note that the exponent of the constant is comparable 
with deg(cr) = k — £ = 2v + s. A factor C £ would not be affordable if £ S> deg(cr); this is 
why the ladder integration had to be done essentially with the precise constant by using 

tirigp . 

Finally, when estimating E(k+ 1), only the last row of M. h = k + 1 G Ii as t is present. 
Since k + 1 is bigger than the largest element of I s U I v , which by Definition 1 10. 31 cannot 
be a covered slope index, we see from the definition of b-^ that b^ = 0, i.e. there is no 
point singularity. Moreover, we have already integrated out \I(\ + \I S \ + 2\I V \ = £ + s + 2v 
variables. Since p = v and I = Ip U I s U I v U I p U Ii as t is a partition, we see that at the last 
step there is only one integration variable left, namely Pk+i- Considering that the (fc+1)- 
th row of M contains only one nontrivial entry and thus [M p + u]fc+i = Pk+\ + 1 , 

the integral from (|10.14j) is simplified to 

E(k + 1) = Fl(/ (fc+1) ,M (fc+1) ,b (A:+1) ) (10.35) 


= \ 2k sup 


Y 


dad/3 


—Y 


1 


dp(pk+i 


1 


< C\ 2k ( d \ log 7/1" 


a - u(p k+ i) - irj\ \(3 - u>(p k +i + u) + ir)\ 


We now combine the inequalities (110.151) . (110.341) and (110.351) . we arrive at 


E(M) < C v+S+1 7 ]-^(^ v +s+\Ius\) ^d(2);+s+ 1 ) | log ^1 2 


where we recall that the general constant C may depend on ipQ. Using the choice of the 
parameters from (I7.25|) and the fact that deg(cr) = k — t = 2v + s > 2, we obtain a total 
A-power 

2(2v + s) — (l + —) (3u + s + |/ us |) — (k + 35)d(2v + s + 1) 

> v + |/ cs | - k( 1 + ^d)( 2v + s) - (2v + s)0(5) 

where we used Lemma 110.41 and v > 1 to estimate ^(3u + s + |/ us |) < 2v + s and 
2v + s + 1 < |(2u + s). Using Lemma 110.41 again, we see that v + |/ cs | > |(2u + s ), i.e. 

E(M) < C ^-( 1+ % d ) K ~ 0{S) y CS(a) \logri\ 2 

if the exponent is positive, i.e. if k < q 2 9c1 and 6 < 5(k) is sufficiently small. This 
completes the proof of Lemma 110.11 □ 


63 


























A Proof of Lemma 110.5 . 

We can replace uj{p) with e(p) in (110.271) by using a straightforward resolvent expansion: 


1 


< 


1 


| a — Lo(p) + ip\ | a — A 2 0(a) — e(p) + ip \ 


1 + 


CX 2 \a — e(p)\ 1 / 2 
| a — Lo(p) + ip | 


< 


C 


|o: — e(p) + ip | 

(A.l) 

with a = a — A 2 Re 0(a). We used the boundedness and the Holder continuity of 0 
(13.51) . Therefore the proof of Lemma 110.51 is reduced to 


Lemma A.l For any |g| < C A 1 

h ■= f 


d nip) 


h := 


a - e{p) + ip | |(3 - e(p + q) + ip\ 
d fi{p) 


< 


C( d 3 | logr/| 2 


(A.2) 


|a - eijp) + ip\ |/9 — e(p + q) + ip\ |||p - r 
uniformly in r, a, (3. 


1 < Cp V 2 C d 3 | log ?7| 2 


Proof of Lemma \A.l[ The bound on I\ follows from a direct calculation and |g| < 
CA" 1 


h < 


l 


dc 


c u d ~ l du r 1 _ 

o |a — u 2 /2 + ip\ 12/3 — (u 2 + q 2 ) - 2|(;|uc| + p 


< 


C( d 3 | logp\- 

ki 


(A.4) 


If p < |(?|, then we use |||g||| < |g| to obtain (IA.2D . If |g| < p, then we use Schwarz 
inequality to separate the denominators and use (13.14|1 to conclude the proof of (IA.2I) . 
To prove (IA.3D . we first establish the following bound uniformly in a: 


J : = 


d Pi?) 


1 


| a — e(p) + ip | |||p — r|| 
that follows by a direct calculation 


J < C 


L 


** u d 3 du 
o \a-u 2 /2 + ip\ 

,f 2 

d- 3 , /-d-2\ 


1 + 


< CC d 2 \ logr/| 


dc 


(A.5) 


< c{ c + c ) 


l |u 2 + r 2 — 2|?’|uc 
dv 


1/2 


<C( d ~ 2 \logp\ 


Jo |a - v + ip | 

with u = \p\, v = u 2 /2 using it 2 + r 2 — 2|r|uc > |u| 2 |l — c 2 | for |c| < 1. 
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We can assume that |g| > p, otherwise we can estimate the /3-denominator in (1A.3D 
trivially by r/ _1 and we can conclude with (1A.5I) . We then distinguish two regimes. If 
|||p — r||| > ry 1 / 2 , then we estimate |||p — r||| trivially and we use (1A.2D . 

Now let |||p — r||| < p 1 / 2 . We split this regime to two subregimes, |cy| < 2|p| and 
\q\ < 2\p + q\, the union of whose clearly cover all values of p. 

In the regime, where |p| < 2|p|, we estimate the square root of the /3-denominator 
trivially and use a Schwarz inequality to separate the remaining (5 denominator from the 
point singularity. The corresponding contribution can be estimated by 

1/2 f l(|g|<2|p|) 

J \ol - e(p) + ir )| 


|/3 - e(p + q) + ip | 


+ 


1 (|||p — r I < p 1 / 2 ) 

Ip - r 1 2 


dp(p) 


The first term was already estimated in (1A.2I) . The second term is bounded by the 
co-area formula by 


Cp 


- 1/2 


L 


J n da 


(\ q \/ 2) 2 |a — a + ip\ \a\ 1 / 2 


with := 


1 (|||p — r 11 < p 1//2 ) 

lib - r I 12 


p) , 


where S a := {p : e(p) = a} and dz^(p) being the surface measure. Clearly J a < |logp|. 
and we obtain the estimate CV/ -1 / 2 ! logr/| 2 /|g|. 

In the regime where |g| < 2|p + ^1, we shift p —> p + q, r —> r — q and interchange 
the role of the a and /3 denominators in the above proof. This completes the proof of 
Lemma IA.ll □. 
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